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Abstract 



The response of vacuum to the presence of external conditions is the subject of 
this work. We consider a generalization of the Casimir effect in the presence 
of curved boundaries on which a sharp potential is concentrated. The profile 
of the potential is a delta function, which has some features in common with 
a hard boundary and some with a smooth background field. The boundaries 
investigated are: i) a spherical shell, ii) a cylindrical shell, iii) a magnetic flux 
tube. The vacuum energy is calculated by means of the Jost function of the 
scattering problem related to the field equation. The energy is then renormalized 
by means of a zeta functional approach adopting the heat-kernel expansion. The 
heat kernel coefficients are calculated and a discussion of the UV-divergences of 
the model is presented. The renormalized vacuum energy is then numerically 
studied and plotted. The sign of E^^"' is found to be always negative in the case 
of the cylindrical shell, while in the case of the spherical shell and of the magnetic 
flux tube the sign depends on the value of the radius. 
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Chapter 1 
Introduction 



The investigation of vacuum energies in the context of quantum field theory under 
external conditions plays an important role in several areas of modern theoretical 
physics. The external conditions can be given by classical fields, like gravita- 
tional fields |l|, 1^, electromagnetic fields, monopoles0, || and more exotic fields 
like sphalerons |^ and electroweak skyrmions |P. External condition may also 
be given by the boundaries of a quantization volume, they can be surfaces with 
an arbitrary shape on which the zero point fluctuations are forced to obey some 
boundary conditions. In the latter case the calculation of the vacuum energy is 
simple, at least for the geometries for which separation of the variables occurs. 
The most elementary boundary conditions are that of unpenetrable surfaces con- 
fining a quantum field in its vacuum state. In the case of two parallel plates in 
the electromagnetic vacuum the result is the well known Casimir effect, which 
has undergone a number of experimental verifications with increasing precision 
0, H, |l^ . For other fields the unpenetrable surface is represented by Dirichlet 
boundary condition in the scalar field case and bag boundary conditions for the 
spinor field case, which find application in the bag model of quantum Chromo- 
dynamics [pT] , |T^. A natural generalization of the Casimir effect are situations in 
which the boundaries become transparent at higher frequencies. An approach to 
this generalization with "softened" boundaries was proposed in |T^, . This 
type of problems are interesting in order to study the response of vacuum to the 
presence of non-ideal boundaries, but also in view of more general investigations 
in the background of a smooth external potential. The latter situation is the 
most difficult to approach. Usually one can perform analytical calculations only 
for some special examples like that of squared potential profiles. Some work in 
this direction has been done in |jl6| , where the vacuum energy of a scalar field in 



the background of a square well potential and of a piecewise oscillatory poten- 
tial were calculated in the one dimensional case. Some numerical evaluations for 



smooth potentials in three dimensions are also known [0]. Given the difficulty to 
solve analytically the problem, it is desirable to investigate a simple model, which 
allows for an explicit solution, while keeping the essential feature of the trans- 
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parency of the boundary unchanged. In the present work we propose a study of 
the vacuum energy in the presence of singular potentials given by delta functions. 
Delta functions are known to be a good idealization of a strong potential concen- 



trated in a very small region [|T8], 0. We will use it to describe thin walls which 
show transmission properties. Such walls can be called "semi-transparent". In 
the context of Casimir like problems a calculation with delta functions already 



exist however only for the case of parallel semi-transparent plates. In the 



present work we will consider more complicated geometries with spherical and 
cylindrical symmetry. In the third chapter a computation of the vacuum energy 
of a massive scalar field in the background of a semi-transparent spherical shell is 
performed. In the fourth chapter we compute the vacuum energy for a cylindrical 
shell. In the fifth and sixth chapter we extend the analysis to a scalar and to a 
spinor field in the presence of magnetic strings. Each of these situations has its 
own history which we describe briefly in the following. 

After the first calculation of H.B.G. Casimir in 1948 ||2T| showed that two 



parallel mirrors in vacuum attract each oder with a force depending on the fourth 
power of their separation, the interest of the community turned to the problem 
of the conducting spherical shell in the electromagnetic vacuum. It was expected 
that the Casimir force would be attractive, since the spherical surface can ideally 
be divided into an infinite number of plane parallel surfaces. Following this idea, 
Casimir suggested in 1953 that an electron might be envisioned as a spherical 
shell of charge e and that the attractive force associated with the electromagnetic 
vacuum might exactly balance the outward repulsive force associated with the 
electrostatic self-energy. The purpose was to explain by means of the zero point 
fluctuations the stability of a charged particle. This fascinating model failed in 
1968, when Boyer showed that the Casimir stress on a perfectly unpenetrable 
spherical surface is positive and that the vacuum tends to expand the shell. The 



result of Boyer, confirmed later by other authors [22, E3l El], put for the first 



time under the attention the strong dependence of the vacuum energy on the 
geometry of the system. Moreover, the work of Boyer showed that the simple 
subtraction of the empty Minkowski space is in general not sufficient to yield a 
finite result. New divergences, associated with the curved surfaces, appear in the 
calculations. The treatment of this divergences needed new and more efficient 



techniques which have been developed only in the last decades [^, ^T\, ^8 
and which are still susceptible of improvements. 

Unfortunately we have no general rule to predict the sign of the Casimir 
energy for a certain boundary shape. The research has gone forward in the last 
50 years with investigations of special cases, from which we try to extract as more 
physical information as possible. A deep understanding of the phenomenology of 
quantum vacuum perturbed by external conditions is still far from being reached. 

Among the special cases, which allow for an analytical calculation, is the cylin- 
drical configuration. Since the separation of the variable is possible in the wave 
equation with cylindrical coordinates and the solution for the field is known to be 
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given by Bessel functions, the summation over the modes can be performed. The 
first complete calculation, realized in 1978 by DeRaad and Milton, ^9) showed 
that the energy of a conducting cylindrical shell is negative. While the inves- 
tigation has turned recently to the problem of compact dielectric cylinders and 



spheres [0, |3T], a treatment of the semi-hard hollow cylinder and sphere 
is still missing. Chapter 3 and 4 of our work, are devoted to the calculations for 
this kind of shells. 

There is an evident analogy between a cylindrical material boundary and a 
magnetic flux tube. A study of the ground state energy in the background of a 
magnetic string is motivated not only by the necessity of completing the class of 
problems mentioned above (for instance with a multiple comparison of the sign 
of the energy), but also in order to gain some insight in to the "interaction" of 
vacuum with purely magnetic fields. Spinor fields in the background of homoge- 



neous magnetic strings are a well know in the literature |3^, |^. Computations 
of the ground state energy were performed by several authors [|36|, ^ which 
extensively considered the case of a string with infinitesimal radius. This model 
has the difficulty that the energy density per unit volume cannot be integrated to 
give the energy density per unit length of the string. Another matter of debate 
is the stability of the string in vacuum: if there exist a value of the radius for 
which the vacuum energy is minimized. Calculation for an homogeneous mag- 



netic field inside a flux tube p9[ and for a planar field with a special profile in 



the x-direction and homogeneous in the y-direction |^ showed the ground state 
energy to be without minimum. In order to shed some light into these topics, 
we propose here a calculation of the vacuum energy of a magnetic flux tube with 
finite radius and with an inhomogeneous potential given by a delta function. The 
calculation is performed for a massive scalar field (chapter 5) and for a massive 
spinor field (chapter 6). 

The technique we adopt for the computation of the vacuum energy makes use 
of the Jost function of the scattering problem related to the background. The 
ground state energy 



■ (n) 



where the uj(^n) are the discrete values of the Hamilton operator associated with 
the background B, can be transformed into an integral over the momentum k 
by means of the Cauchy theorem and with the deformation of the integration 
contour on the imaginary axis. Then we can writeQ 

Eb = -S^^^sV- rdk{e-m'y/'-^^\nm) , (1.2) 



"'^This equation is valid in a spherically symmetric background. For a cylindrical geometry 
the equation changes slightly. 
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here J dk and Yli ^^re respectively the integration over the momentum k and 
the summation over all the remaining quantum numbers /; m is the mass of the 
quantum field, s is a regularization parameter and fi a mass parameter which one 
introduces to maintain the correct dimensions of the energy. The Jost function 
fi{ik) is given for each problem by the coefficients of the regular solution of 
the wave equation. A nice feature of this integral representation is that the 
contribution of the possible bound states is "automatically" included. 

The renormalization of (1.2) is performed with the zeta functional method (for 
a review see In this approach the divergences are completely described by 

the lowest heat-kernel coefficients associated with the field equation. We expand 
the regularized vacuum energy in powers of the mass of the field. This expansion 
contains the heat-kernel coefficients. The renormalization is then performed by 
subtraction of the pole terms and of the terms proportional to the non-negative 
powers of the mass. The reason for this, is that the renormalized vacuum energy 
must satisfy the condition 

= for m^oo (1.3) 

which is natural under the physical point of view. This condition defines the 
renormalized vacuum energy removing the ambiguity introduced with the mass 
parameter /x. 

For each of the configurations considered in this work we will calculate the 
relevant heat-kernel coefficients. Their importance resides not only in the calcu- 
lation of the UV-divergences, but also in the study of the asymptotic behaviour 
of the energy, which is provided by the lowest half-integer heat-kernel coefficients. 

Other mathematical tools are employed in the course of the calculations, like 
the introduction of the uniform asymptotic expansion of the Jost function to 
perform the analytical continuation s — > 0; and the application of the Abel- 



Plana formula for half- integer variables The latter technique allows for an 



analytical elaboration of the vacuum energy toward a more compact form. After 
the renormalization the vacuum energy is evaluated numerically as a function of 
the relevant parameters of the background, namely the radius of the shell and 
the coupling. At the end of each chapter we display plots of the renormalized 
vacuum energy. The plots are briefiy discussed in a special section for each of the 
4 problems. The global results of the work are then summarized and commented 
in the conclusion. We would like to begin with a short introductory description 
of vacuum in QFT and of the zeta functional renormalization. 
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Chapter 2 

Vacuum energy and 
renormalization 



2.1 Vacuum energy in quantum field theory 

Quantum field theory is a formulation of the laws of nature which allows to predict 
and describe the phenomena of particle physics. The canonical quantization of a 
field (or second quantization) follows a standard procedure, which we summarize 
briefly in the following for the case of a massive scalar field. The Lagrangian 
density of a real field (f){x) with spin and mass m, is 

where natural units h = c = 1 are used. From the Euler- Lagrange equation 
1^ = -^Q^§r^ one arrives at the Klein-Gordon equation 



,2\ 



(□ + m')(/)(x) = 0. (2.2) 
With the definition of the canonical conjugate field 

the Hamiltonian density in terms of canonical variables reads 

n = 7r(x)0-£ = ^(7r'(^) + (V(/)(x))2 + mV(^)) ■ (2-4) 

In order to quantize the system the canonical variables are transformed into the 
operators and tt obeying the commutation relations at equal time 

[^{x,t),TT{x',t)] = iS^^\x-x') (2.5) 
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= [n{x,t),7r{x',t)] = 0. (2.6) 
Then, the quantized Hamiltonian reads 

H = j d^x]^(T:'^{x,t) + (V0(f,t))2 + m202(f,t)) (2.7) 

Now, we decompose the field and its canonical conjugate using a basis of eigen- 
functions e*^^ and introducing the operators dp and at 

= / ^== («pe'^^ + aje-^n (2.8) 

rr(f, t) = j i-tu,) {d,e^^^ - aje-^^') , (2.9) 

where Up = \/^^\-~m? is the energy of a particle with momentum p and the 
operators dp and at obey the commutation relations 



[dp, dp, 



S(^){p-p) (2.10) 



[dp, dp] = [at,aj^] = , (2.11) 

which can be proven by means of ( p.6|) and (|2.5|) . The relations ( p.lO|) and 
( |2.11| ) represent a Fock space in which the dp and at play the role of creation 
and annihilation operators respectively. Inserting ( p.8|) and (|2.9| ) in { \2.7\) we can 
express the Hamiltonian in terms of these operators 

H = ^ j d^ptupididp + Op-at) . (2.12) 

To give a physical meaning to this expression we enclose the system in a large 
quantization box of dimensions LxLxL on whose surface the field obeys periodic 
boundary conditions. Then the integration over the momentum in the continuous 
is replaced by a sum over the spectrum of the discrete eigenvalues or "modes" of 
the field and the Hamiltonian becomes 

H = ^^u;j(ajaj + a^aj) . (2.13) 
i 

The commutation relation ( |2.10| ) becomes 

[d„d],] = (2.14) 
With the help of (|2.14| ) the Hamiltonian is rewritten as 

H = J2uj,id]d, + ^) . (2.15) 

j 

8 



In this form the Hamilton operator can be easily interpreted: in each sate j there 
are some particles whose number is given by fij = OjCij and the energy quantum 
is ujj . Moreover it appears in ( |2.15| ) a contribution ^Uj which is independent of 



the occupation number. It is the so called zero point energy or vacuum energy 
or, if one wishes, ground state energy, in analogy with the lowest energy level 
of an harmonic oscillator in quantum mechanics. This contribution is not a 
peculiarity of the Klein-Gordon field. It appears in the Hamilton operator of each 
quantum field (complex scalar fields, Dirac fields. Maxwell fields, gauge fields in 
Chromodynamics and so on). One must only be careful with the commutation 



relation ( 2.14 ), which changes in each case. For instance in the quantization of 



a spinor field the addend | in ( p.l5|) appears with a sign minus, because of the 



use of anticommutators in the Fermi-Dirac statistics. Furthermore the index j 
must take into account all oder possible degree of freedom of the field. So it is 
usually replaced by a subscript (n) which includes all the quantum numbers in a 
compact form. The vacuum energy 



En 



lY.^in) (2.16) 



2 

in) 



is a divergent quantity. Since physically we can only observe energy differences, 
the quantity ( |2.16| ) as a constant of the energy, cannot be physically observed. 
However a variation of Eq caused by some external conditions can be measured in 
laboratory. We must not forget that result (|2.15|) is obtained by imposing periodic 



boundary conditions on the surface of a quantization box of size L, which is a 
fictitious boundary. If we introduce a physical boundary on which the field is 
forced to obey some conditions, the modes of the field will change and with them 
the vacuum energy. The most eminent example of this situation is the Casimir 
effect. This is the attraction of two conducting parallel plates in the vacuum of 
the electromagnetic field. The force is caused by a difference in the spectrum 
of the electromagnetic wave vector k = {k^, ky, k^) inside and outside the plates. 
Because of the boundary conditions: E\\ = 0, B± = for an electromagnetic 
wave near a perfect conductor, the component of k perpendicular to the plates 
can only take the discrete values: kz = nn/d inside the cavity, where d is the 
distance of the plates and n = 1,2, ... All other components inside and outside 
the cavity have a number of allowed values which is so large that can be well 
approximated by a continuous spectrum. Then we have 

AEo = Eip - 



j dk.j.dky\jkl + k. 



d'kJkl + kl + kl , (2.17) 
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where the two independent polarisation states have been taken into account. This 
difference was computed by Casimir |^ with the help of the Euler-McLaurin 



formula (see [^], equation 3.6.28). The result is a negative energy: AEq = 
~7^cP 'which gives rise to an attractive force between the plates. 

The conducting plates are just an example of a wide class of boundaries which 
influence the vacuum. For a scalar field the boundary conditions can be Dirichlet 
boundary conditions 

Hx)U=.B = 0, (2.18) 

where a;^ is the value of the spatial coordinates at the surface of the boundary, 
or Robin boundary conditions 

dn{u{x)<j){x))\,=,^ = , (2.19) 

where 9„ is the derivative in the direction normal to the surface of the boundary 
and u{x) is an arbitrary function of the coordinates. In the case u{x) = 1 we have 
the Neuman boundary conditions. For a spinor field local boundary conditions 
are usually studied 

(l + n7)V^(x)U=,.5 = , (2.20) 

like in the case of the bag-model for the hadrons. The variation of the spectrum 
(|2.16|) can also be generated by the presence of a classical background field. This 
is usually introduced with a potential V{x) which can couple to the mass in the 
case of a Newtonian background field or can enter the covariant derivative, as in 
the case of a gauge potential A^{x) which yields the operator ^ + iA^{x) 
in the field equation. In problems with scalar fields the potential can be given in 
the simple form a/R where a is a generic coupling constant and i? is a geometric 
parameter. In the context of general relativity the metric Qf^i, can also play the 
role of the background potential. 

In each of the above mentioned situations the calculation of the zero point 
energy requires the knowledge of the spectrum ujB{n) where the label B represents 
here a boundary condition or a background potential. The calculation of the ujB(n) 
can be a hard task, when the geometry is non-planar and when the potential is a 
smooth function. The other difficulty related with the calculation of the vacuum 
energy, is the renormalization. In the Casimir effect, the subtraction performed in 
( |2.17D is just an elementary example of renormalization. Namely it is subtracted 
the energy of the empty unbounded space (Minkowski space contribution). This 
simple procedure does not give a finite result in spherical or cylindrical config- 
urations. For these situations the isolation of the UV-divergences is performed 
by asymptotic expansions of the vacuum energy in terms of a cut-off function. 
In recent years the zeta functional technique is used for the renormalization. 
We briefly describe the technique in the next section. For a more exhaustive 
treatment the reader is referred to ^ . 
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2.2 The zeta functional method 



Zeta functions provide an elegant re- formulation of the sum (|2.16| ). The sum is 



regularized by a complex parameter s which must vanish after the renormaliza- 
tion. The regularized energy has the form of a zeta function, whose poles are well 
known and can be subtracted to yield a finite energy. Moreover, a zeta function 
has a variety of interesting properties, like the integral representation and the 
heat-kernel expansion, which allow for easy mode summation and for a deeper 
insight into the structure of the UV-divergences of the investigated model. 

2.2.1 Integral representation 

Let us start with a general review of the analytical properties of zeta functions. 
The most simple zeta function is 

C(^) = E;^- (2.21) 

n 

This is commonly known as the Riemann zeta function Cr{s), since Riemann 
studied the general case of a complex variable s and found the analytical con- 
tinuation in the region < 1, where representation ( |2.21| ) does not converge. 



The analytical continuation is of interest also for physical applications. In quan- 
tum field theory one generally adopts the zeta function related to a differential 
operator. If one has to consider a field cf) and the wave equation 

D4>n = en^n , (2.22) 

where D is a generic differential operator (like a Laplace operator) and e„ are its 
eigenvalues, then the zeta function of the operator D with spectrum T> is defined 
as g 

Cnis) = 5^(6„)-^ . (2.23) 

The Hamilton operator of a system is then given by the zeta function of the 
corresponding differential operator with a simple function of s in the argument. 
For instance, in the problem of a spherical shell which we will analyze in the 
next chapter, the regularized vacuum energy is given by ('(s — 1/2), while in the 
case of the cylindrical shell the relevant zeta function will be ({s — 1). They are 
elementary generalizations of the Riemann zeta function. The latter can undergo 
at least two types of transformation into an integral. The first makes use of the 
integral representation of the Gamma function 

POO 

T{s) = dt f-^e-' , (2.24) 
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from which it follows 



n 



and 



oo „ 
n=l 



Tis) 



CO +S-1 



s~l „—nt 



dt f~'e 



dt 



-nt 



r s) 



oo I 

dt -=-rT 

r(s) e*-l 



(2.25) 



(2.26) 



which in the case of the function Cd{s) becomes 



Cd{s) 



oo 

dt — - 

r(.) 



(2.27) 



The second representation makes use of the Cauchy theorem, which transforms 
the sum ( p.21| ) into a contour integral 



dt- 



,2nit 



- 1 



(2.28) 



The contour 7 encloses, in the complex plane, all the positive integers lying on 
the axis 3fJt, as it is shown in the figure below. Each integer is a pole of the 
integrand in ( |2.28D and the residue contributes to the initial sum ( |2.21| ). 



"^t 
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Figure 2.1: The integration contour of (2. 



The integral 



is convergent only for 3fJs > 1. We are interested in a 



representation in which the continuation for small s is possible. To this aim we 
split the integration contour in two pieces: an upper contour 7_|_ for Qt > and a 
lower contour 7. for "^t < 0. With the help of: l/(e2*''*-l) = -1 + 1/(6-^^'^* - 1), 
we get 



dt r 



4 



l-e 



+ 



dt t~ 



1 



(2.29) 



The integration over 7+ is split into two parts, the first of which can be explicitly 
calculated and it gives e^~*/(l— s), where e G (0, 1) is the point at which the initial 
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contour 7 crosses the real axis. In the second part we deform the integration path 
until it runs along to the positive imaginary axis from e to ioo. On this path 
the integrand function falls exponentially. In the integral over 7„ we deform the 
contour downwards until it runs from e to —ioo. By this means we find a result 
in which the analytical continuation to s < 1 is possible. After setting e = we 
find 

/•oo J.— s 

Cr{s) = 2sin(s7r/2) / dt^^^ . (2.30) 
Jo e — i 

This expression can be generalized for a function Cd{s), as we will see in each of 
the following chapters, providing a useful tool for the calculation of the ground 
state energy. 



2.2.2 Renormalization via heat-kernel expansion 

The regularized vacuum energy must be at the end a finite quantity for s = 0. 
To achieve this result, the divergent contribution to the zeta function must be 
isolated. This is done by means of the so called global heat-kernel K{t), given by 

00 

K{t) ^ J]e-"*. (2.31) 

n=l 



This expression appears already in ( |2.26 ). It its expansion for small t which is of 
interest : 

00 

Here d is the dimension of the manifold and the index j assumes also fractionary 
values in the case of a manifold with boundary. The coefficients aj are called the 
heat-kernel coefficients. If we insert ( |2.32| ) in ( p^.26 ), we get a formula in which the 



divergences appear associated with the lowest heat-kernel coefficients. To better 
enlighten the problem and also to restrict us to the renormalization procedure 
which we will adopt in this work, let us consider a physical system given by a 
massiveQ scalar field and by a classical background described by a potential Vh{x) 
concentrated on a surface with arbitrary shape. The ground state energy of the 
system is given by 

= \T. \/pU + ' (2.33) 

in) 



where ypfn) + the energy eigenvalues of a particle of the field and m is 

the mass of the field. We regularize this expression with the introduction of a 
"'^We will work throughout this thesis with massive fields. 



13 



parameter s, to be put to zero after the renormalization, and we rewrite ( |2.33| ) 
in terms of a zeta function 



= lf^''J2^pU+^'y^"" = lf^''(b{s-l/2) , (2.34) 

(n) 

where fi is an arbitrary parameter with the dimension of a mass. Then we have 

^^^'X^{s - 1/2) = lf^''Y^^^Y/2) J ^'"'^'^''(^) • (2.35) 
With the help of 

(n) i 

we can expand the zeta function in a series in which the poles 1/s, contributing 
to the divergences, are well isolated 

^^ = ^^^(^-^/2) - -64^U + ^^7^-2j^°-24^"V. 

/ 1 , 4u^ ^ \ m 

+32^ u+^"^" 'r^ +T^"^/^ 



1 fi 



327r2 



fi + ln^- 2^2 +0{m-^). (2.37) 



Here the terms of order 0{s) have been dropped. The renormalization of the 
ground state energy is then performed by simple subtraction of the pole terms 
~ 1/s, corresponding to the lowest heat-kernel coefficients. This subtraction is 
physically interpreted as a re-definition of the classical parameters of the system. 
In fact the heat-kernel coefficients depend on the geometric features of back- 
ground. For instance, if we calculate the vacuum energy inside and outside a 
spherical boundary, the coefficients aj will depend on powers of the radius R of 
the sphere. Now, it is possible to express also the classical energy of the sphere 
in powers of the radius 

Ef'^'' = (^pV + aS + FR + k + ^^ , (2.38) 

where V is the volume of the sphere, p is the pressure, a is the surface tension 
and F, k and h are other parameters with no special names. Then, the renor- 
malization of the vacuum energy requires a renormalization of the parameters 



^Such a definition appeared for the first time in P7[, see also section 3.1 for other details. 
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defined in (UM 



9 



327?^ V s m2 



h/R - Vi?-^(^ + ln^- 2) a,. (2.39) 



By means of tliis procedure, the total energy of tlie system remains unchanged. 
The transition from ( 2.38| ) to ( p.39 ) does not correspond to a physical phe- 



nomenon. The parameters ( ^.381 ) are "naked" quantities which do not take in to 
account quantum corrections. In an experimental measurement only the parame- 
ters given in (|2.39|) are observed. An important question is if the renormalization 



procedure described here is unique. In the case of a massive field, it is not clear 
which value the mass parameter fi (introduced for dimensional reasons) should 
take. Generally this ambiguity is removed by the requirement that the renor- 
malized vacuum energy must only depend on negative powers of the mass, which 
is sufficient to fix univocally the meaning of the vacuum energy. In this case 
the terms proportional to and in (|2.37| ) must also be subtracted. They 



correspond to the heat-kernel coefficients with fractionary index. 

The heat kernel coefficients play a crucial role in this technique. An explicit 
knowledge of the lowest coefficients fully solves the renormalization problem for a 
given configuration. The coefficients ai, 02 and 03 are known, at least in principle, 
for a generic background potential Vh- 

Oo = / d^x , 



ci = — / d^xVbix 



a2 = d^xVi^ix) . (2.40) 

The first coefficient represents the contribution of the flat empty space, which, 
as we saw, is subtracted for the derivation of the Casimir effect. In the simplest 
example of two parallel plates all other coefficients vanish. In general the coef- 
ficients Oi and 02 are not zero and must be calculated. Unfortunately formula 
( |2.40D cannot always be used. For instance, in the problems that we analyze in 



this thesis the background potential contains a delta function and the equation 
for 02 would be ill-defined. One last comment must be made about the limit of a 
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massless quantum field. In this case equation ( p.37| ) does not apply and another 
type of expansion is needed. Furthermore the requirement that the vacuum en- 
ergy vanishes for m — > cxo would fail and the renormalized energy would not have 
a unique meaning. However the present work is restricted to fields with finite 
mass. In the following chapters the mathematical tools summarized here will be 
applied to four different exercises. The vacuum energy will be calculated by mode 
summation with integral representations of the type (|2.30| ). The renormalization 
will follow formula ( |2.37|) . Modifications in the procedure, due to the specific 
features of each exercise, will be displayed and explained in each chapter. 
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Chapter 3 



Vacuum energy for a 
semi-transparent spherical shell 



The interest in the zero point energy for spherically symmetric configurations 
has been renewed in the last years by the phenomenon of sonoluminessenceQ 



which, in the suggestion of some authors could be a macroscopic 

manifestation of vacuum. This hypothesis, which has been lively discussed by 
the community, could not be completely accepted or ruled out, since an exact 
model for the dielectric ball in vacuum is still missing]^. The model of a hollow 
spherical surface has recently attracted the interest of physicists in the context 
of the MIT bag-model for the hadronic structure. In this model, confinement is 
provided by special boundary conditions requiring the vanishing of the quark and 



gluon currents through the boundary (see eq. ( |2.2(J| ), and |n|]). Massive scalar 



fields bounded by a Dirichlet spherical shell were considered recently in ||4^, with 
the result that the sign of the ground state energy is positive. 

In the present chapter we calculate the ground state energy of a massive scalar 
field in the background of spherical shell with delta function. The technique 



developed in for generic background fields is used and explained with details. 
The calculation of the ground state energy is based on the knowledge of the 
Jost function of the associated scattering problem, which is given by a simple 
formula in terms of modified Bessel functions. The renormalization is performed 
via zeta function. The heat-kernel coefficients whose knowledge is required for 



the renormalization, were first derived in Here they are rederived in the 

course of calculations. 



^The emission of short intense pulses of light by collapsing bubbles of air in water [ p4[ . 
^The compact dielectric ball was also in the center of interest at the workshop "Quantum 
field theory under external conditions" held in Leipzig on 14-18 September 1998. 
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3.1 The model 



We want to study the ground state energy (GSE) of the scalar field '^(t,x) to be 
quantized in the background of a potential V{r) concentrated on the surface of 
a sphere. We start with the following field equation^ 

{n + m^ + V{r)) <^{t,x) = 0, (3.1) 



where m is the mass of the field and r = a/x^ + + z^. This equation will 
be rewritten later in polar coordinates. The sphere is a geometrical object with 
radius R and a surface S, to whom it can be associated a classical energy in 
terms of classical parameters as we mentioned in the previous chapter. The total 
energy of the system reads 



Etot — Edass + 



quant 



(^pV + aS + FR + k + ^^ + ' (3.2) 



The classical part of the energy is expressed in a general form in which the depen- 
dence on powers of R is explicit. This definition is suitable for its renormalization, 
it has been introduced in and used in many works concerning the bag model 
and the Casimir energy for fermionic and scalar fields with spherical boundaries 
pUj . The quantum contribution in ( |3.2D is the expression for the vacuum energy 
of a scalar field whose energy eigenvalues are Un- To render the eigenvalues of 
the energy discrete we take temporarily a finite quantization volume with radius 
L > 1. 

The classical shell is static and spherically symmetric. It is described by a 
potential 

V{r) = ^6{r-R), (3.3) 

where a is the strength of the potential. The continuity of the field on the 
boundary will be discussed later. The potential could be expressed also in other 
forms involving the mass, for instance as 

V{r) = a m 6{r — R) 

since both m and R are dimensional parameters. However the choice of R is the 
most natural since the mass concerns the quantum field while the radius concerns 
the background potential which is independent from the field. 

The quantum contribution to the total energy is divergent, for the regular- 
ization we adopt the zeta function technique . We define the regularized ground 
state energy 



E^ 



^E(^(™)+-^')'^">"' (3-4) 



(n) 



•^Natural units are used 



18 



where /x is an arbitrary mass parameter, s is the regularization parameter which 
we will put to zero after renormalization and A(„) are the eigenvalues of the wave 
equation 

[-A + Vir)]^(^n){x) = \ln)nn){x). (3.5) 

Now we introduce a zeta function. The zeta function of the wave operator with 
potential V{r) as defined in ( p. 5] ) is 

Cv{s) = J2(^U + mT- (3.6) 

in) 

We can express the ground state energy in terms of this zeta function 

K = lCv{s-l)f^''. (3.7) 



Using 



= —— / dt f -^e"^*, (3.^ 

r s 



we can write eq.(2.6) in the following form 

1 /"OO 

Cv{s) = ^J dt f-'e-"^'' Yl (3.9) 



that is 



in) 



Cv{s) = dt t^-'e~"^''K{t), (3.10) 

r(s) Jo 

where function K{t) is the heat kernel. Taking its the asymptotic expansion for 
t 0, 

/ . N 3/2 oo 

ir(t) = 5^exp(-Aj„)t)*^° — j = 0,-,l,..., (3.11) 

(n) V ^ / j=o 

and making the substitution s— s^s — l/2in eq. (|3.10|) , we get an expansion of 
the type ( |2.37[ ), in which it is easy to recognize the pole terms. This makes it 
possible to define the total divergent contribution to the ground state energy by 

Tj^div (I 1 \ 

/ 1 , 4u^ \ , m . 



J— (^1 + In 1^ - 2^ ^2. (3.12) 
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The quantities Aj are the heat kernel coefficients of the spherical 5-shell. In 
the definition (|3.12| ) we have included the terms A1/2 and with half integer 



index, which do not contain poles, to satisfy the normalization condition which 
we formulate explicitly below. We define the renormalized zero point energy by 

E;''' = E^-Ei'^. (3.13) 

To keep the total energy of the system unchanged one must add the subtracted 
object i?^™ to the classical energy. Then one has also a definition of a new 
classical energy 

^dass = Eclass + E^^^ . (3-14) 

The transition from E^ass to edass consists in the renormalization of the classical 
parameters in a way described in ( |2.39| ). However in our particular case both the 
Aq and Ai/2 coefficients, corresponding respectively to p and a, will turn out to 
be zeroQ, then only the last three terms in ( ^.39| ) will undergo renormalization. 
Now we have 

EtOT = ^class + -E'^'^" • 



The old classical energy Edass as defined in (|3.2| ) is an unphysical quantity, since 
experimentally we can observe only an energy which includes the vacuum fiuc- 
tuations. The term h/R in ( |2.39|) deserves a particular attention. In fact, in 
the case of a massless quantum field the vacuum energy takes the form ~ 
Therefore the classical and the quantum contributions would not be distinguish- 
able and the calculation of E^^"' would lose its predictive power. This difficulty 
makes it impossible to apply our procedure in the limit m —>■ 0. Furthermore, we 
must note that the ground state energy proposed in ( |3.13|) has not yet a unique 
meaning. For the uniqueness of E^'^"' we impose the normalization condition 

lim E;"" = 0, (3.15) 

which physically means that for a field of infinite mass we have no quantum 
fiuctuations. We fulfil this requirement by subtracting all the contributions in 
£jdiv proportional to non negative powers of the mass. That is we subtract also 
terms with fractionary indexes up to and including the term resulting from the 
heat kernel coefficient A2. The remaining part, containing only negative powers 
of m, will go to zero for m 00. Note that condition (p.l5 ) does not apply to a 
massless field. 



''More exactly, the contribution of the Minkowski space ^0 does not depend on the back- 
ground and can be simply ignored. 
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3.2 Representation of the ground state energy 
in terms of the Jost function 



We adopt the approach appeared for the first time in for the calculation of 
the GSE in the background of a smooth potential. The method for the calculation 
of the heat-kernel coefficients for different boundary conditions was developed in 



an earlier work 53 



With the ansatz of the separation of the variables in spherical coordinates, 
one finds the following radial Schrodinger equation 

1(1 + 1) . X ,9 \ /X / X 



where / is the angular momentum. In the general scattering theory with a con- 



tinuous spectrum p we have the "regular solution" [54| defined as 



where jiipr) is the Riccati Bessel function. The asymptotics of the regular solu- 
tion is expressed in terms of the Jost function fiij)) 

^vAr) ™ I {fi{p)kipr) - fiipfhtipr)) , (3.18) 

where hf{pr) are the Riccati-Hankel functions. Now we examine the field at the 
boundary of our quantization volume. As the potential has a compact support, at 
the boundary, expression ( |3.18| ) becomes an exact equation. It can be considered 
as an equation for the eigenvalues p = Xn,i- Now, taking for instance Dirichlet 
boundary conditions at L: ipp^i{L) = 0, we get 

mkipL) - ft{p)hi{pL)) = 0. (3.19) 

Since eq.( p.l9| ) is satisfied for p = Xi^n, we can rewrite the sum in (|3.4|) as a 
contour integral using the Cauchy theorem 



= /.^^f:(/+i/2) [^{p 



d_ 
dp 



In [fiip)KipL) - ft{p)ht{pL)) , (3.20) 



where the contour 7 encloses all the solutions of eq.( |3.19 ) on the positive real 



p axis and also the bound state solutions in the limit L ^ 00, which lie on 
the imaginary axis. We further simplify eq.( |3.20|) by separating the contour into 



two pieces 71 and 72 and expanding the Hankel functions for large L. Then it 
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is possible to recognize in the integrand a term ipL which corresponds to the 
Minkowski space contribution. This term can be dropped. Now we shift the two 
contours 71 and 72 to the imaginary axis and substitute p —>■ ik. We then obtain 

00 /»CXD O 

^^ = -^/i^^5^(/ + l/2) / dk[e-m']'/'-^^\nMtk). (3.21) 
^ 1=0 ^"^ 

Since at the end our quantization volume will go to infinity { L ^ 00) this equa- 
tion will be independent from the boundary condition chosen for the quantization 
volume. Eq.( p.2lD is a very general and useful representation of the ground state 
energy, where all the information about the background potential is contained in 
fi{ik), possible bound states as well. 

In order to perform the analytical continuation to s = and the subtraction 
proposed in (|3.13| ) we split i?^^" into two suitable parts, which are convergent 
in the limit s — > 0. We obtain this by adding and subtracting the uniform 
asymptotic expansion of the Jost unction (for more details on this procedure see 
ll). We define 

E:;^ = Ef + Eas, (3.22) 
E, = -^/.^^J](/ + i) |^dk[e-m'Y/'-^^[lnM^k)-lnfr{^k)] (3.23) 

and 

E.. = -^/.2^^(/ + l) l^dk[e-m'f'-^^\nfrm-Ei-, (3.24) 

where fi'^{ik) is the uniform asymptotic expansion of the Jost function (taken 
for u and k equally large) which we will take up to the third order in z/ = Z + 1/2. 
In fact, three orders are sufficient to make (|3.23| ) converge and they allow to put 



s = under the sign of the sum and the integral. Higher orders could be included 
to speed up the convergence, however, the final result E^^"' remains unchanged 
in whatever order (> 3) In fj^^{ik) is taken. This is obvious, since the quantity 
subtracted in the integrand of (|3.23| ) is added again in (|3.22|) . 



Now we need the Jost function corresponding to our scattering problem, so 
we turn to study the field at the surface of the shell. 

3.3 Jost function of the spherical (5-shell 

The initial field equation 



l{l + 1) 



dr'^ 



- + AlJ ^n,K^) = 0, (3.25) 
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valid for — oo < r < oo can be divided into two parts: an equation for the free 
field 

at r^R — > (n + m^) (^(t,f) = 0, (3.26) 

and an equation for the field on the shell, which includes the matching conditions 
based on the continuity of the field at the boundary 

_ J continuous. , . 

= ^ \ ip'{R + 0) - ip'{R - 0) = '> 

where the prime indicates derivative with respect to r. We take the regular 
solution for the delta potential 

^kAr) = Ji{kR)Q{R -r)+'- [fiik)h-ikR) - ft{k)hi{kR)) 0(r - R) (3.28) 

consisting of two pieces inside and outside the radius i?, respectively. As above, 
ji{kR) is the Riccati-Bessel function and hf{kR) are the Riccati-Hankel functions. 
Combining eq. (|3.27|) with eq. (|3.28|) we get 



3i{kR) = f ^mhiikR) - ftmrm^ . 

I UkR) = k U (Mk)h'r{kR) - f*{k)h[+ikR)) -jlikR)'' 



We solve for fi{k), keeping in mind that the Wronskian determinant of hf is 2i. 
We find 

m = i (-2z(-l) + 2i^3i{kR)ht{kR)) (3.30) 

or 

rv 

fi{k) = l + —ji{kR)hi{kR). (3.31) 

For the Jost function on the imaginary axis we get 

Uiik) = 1 + ah{kR)Ky{kR), (3.32) 

which is in terms of the modified Bessel functions ly and K^, where v = I + 1/2 . 
We need also the asymptotics of the Jost function: f^'^iik), or more exactly the 
logarithm of fl^^{ik). The expansion of the product of the two Bessel functions 
in ( |3.32|) , for k and u equally large, is easily obtained with the help of |Q. Then 



we find the needed asymptotics as a sum of negative powers of u with coefficients 
Xj^n depending on a. We define 

n=l j 

at a 3a 5a f 



with t = l/^l + Now, inserting (jpl) and (lOSI) in (Q) and dpI) the 



renormalized ground state energy E^^"' = Ef + Eas can be calculated. 
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3.4 Analytical simplification and heat-kernel co- 
efficients 



Let us work on Eas- We transform the sum over / into an integral with the help 
of the known Abel-Plana formula 

|.,4,.f...,,H-f^^M^. ,3.34) 



In our case we have 



/oo a 
dk V \e - m2]V'''^ In (3.35) 

which analytically satisfies the validity conditions for eq.( p.34| ). Eas is split into 
two addenda: 

e£) = -^^- rrf.F(.) (3.36) 

and 

71 Jo l + e^""" i 
First we calculate E^J; for the k and i/-integrations we use the formula 

Jo Jm dku^ 2(mi?)"-2r(|) 

(3.38) 

then 

EW = -""^^^'^ duu P dk[k^-mY"~^^\nn\ik)) (3.39) 

^l-2s^2s^^ (^i^)2-n r(| - g)r(i + i^)T{s + ^) 



Here, inserting the coefficients of ( p.33| ) and expanding up to the first order in s 
all the terms which depend on the renormalization parameter we get 

E^as = - + In -h- - 2 + In -h- - 1 H . 

967ri? \s J Svr \s J 16 

The terms poles contribute to the divergence of the energy. They are used to 
calculate the heat-kernel coefficients in ( 3.12| ) and will disappear after the sub- 
traction of -E^™. The term proportional to m corresponds to the ^3/2 term of the 
heat kernel expansion; although this term generates no divergence it will be as 
well subtracted because of the normalization condition. 
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(2) 

Now we calculate EaJ ■ The integration over k is carried out with the formula 



r) 

dk -ml^/'-^^t^' 
ok 



-m 



,_2,r(f-s)r(.+ 



2 ' 



\mR) 



r(|) 



j-1 • 



(3.40) 



We obtain 



TT 



3,n 



r(i; 



[IV 



This expression can be transformed into 

4/i2 



- + 



^ ^ KmR) 



. (3.41) 



+ 



1 

32^2 Vs 
2a^ + 3« 



+ In 



Q! Q. 

2 1 - hi{Rm) - — — 62(i??Ti) 



bs{Rm) 



ex ex 

hi{Rm) + — — h^{Rm) 



(3.42) 



487ri? ' SttR ' ' 48nR 

where integration by parts was used and the following functions containing the 
z^-integrations of eq. (|3.41|) are introduced: 

,2 



bi{x) 

63 (x) 

64 (x) 

65 (x) 



du 



dv- 



\ ln|l-^|; 



1 + e 



1 _ £1 



/°°.w„|i-^l(- 



o27ri/ / ' 



1 (I 



du\n\l-—\ - - 



V \ v \ 1 + e 



(3.43) 



We see that also EaJ has a pole of the form for j = n = 1. This pole contributes 
to the heat kernel coefficient A2. 

Now, since Ej contains no poles, we are able to write down the complete heat 
kernel coefficients Aj, up to the order j < 2: 



A 



1/2 



0: 



2 a 

Ai = -A-nRa; = vr^/^a^; = --tt— 



25 



This coefficients are the same as in paper After performing the subtraction 

E^J cancels completely and only e!^J (without its divergent portion) contributes 
to the total energy. So we have finally, 

Eas \s-Q = oARm) — OoiRm) H b-AHm) 

' " 2'kR ^ ' 8R^m ^ ' 48nR ^ ' 

a a 

hi{Rm) + — — h{Rm) . (3.44) 



3.5 Asymptotics of E, 



'as 

It is interesting to check (analytically as far as possible) the behaviour of Eas for 
small and for large values of R. To this aim, we first calculate the corresponding 
asymptotics of the functions hn{x). 
In the limit — >^ we find: 

YvsxvhiiRm) ~ — ln(i?m) + Ci ; 

\\mh2iRm) ~ R^n? + C2 ; 

lim b^^Rm) ~ ln(i?m) + C3 ; 

R — *0 

\imhi{Rm) ~ 2\n{Rm) + C4 ; 

R — >0 

\mvh{Rm) ~ 81n(i?m) + C5 , (3.45) 

i? — 

where the Cn's are numbers resulting from the ^/-integrations. Then we have 

R^o IGnR 247ri? V Rm J 

where Cq = {-8Ci + C3 - C4 + C5/3) ~ 0.224 and C3 ~ 1.96. Therefore, 
the asymptotic part of the energy behaves logarithmically for small radii of the 
spherical shell. This behaviour will be observed numerically in the complete 
renormalized ground state energy. 



For R ^ 00 we have 



7 1 



lim bi(Rm) ~ ^ ^ 

lim boiRm) ~ — : 

lim b2,{Rm) ~ — „ „„ ; 

7 1 



lim bARm) ~ , , 

i?-oo ^ 480m4i?4 ' 

31 16 

lim 65(i?m) ~ — , (3.47) 

iJ-oo ^ 16128 m6i?6 ' V ; 
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and we find 



lim Eas ~ 



a 



384 m2i?2 ■ 



(3.48) 



3.6 Numerical results 

To numerically study Eas we rewrite it in the form 

1 



E„ 



16nR 



(3.49) 



where the three functions Vn{x) are given by 
vi{x) = 

V2ix) = 



lbi{x) + hix) - 2b^{x) + -h{x) 



X 



VsiX 



^63 (x). 



(3.50) 



As it is clear from (p.49|) , for small values of a Eas behaves like function vi(x). 



For large values of a, Eas behaves like function V3{x). The plots of the three Vn{x) 
functions are shown below in Fig. 3.1. In this as in all the following plots m is set 
to be equal 1. For the complete quantum energy we still need the contribution 




Figure 3.1: The curves of the Vn{x) functions. 



Ef. In the expression (p.23|) , after putting s = 0, we integrate by parts and we 
obtain 



Ef 



1=0 •'^ 



k 



{In U{zk)- In f^'{zk))dk. 



(3.51) 
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Figure 3.2: The curve of R ■ Ef{R) for a strength of the potential equal to 1. 
For R = Q the curve converges to a finite value. 



This quantity cannot be further analytically simplified. Below we show (Fig. 3.2) 
a plot oi R- Ef as function of i? for a = 1. 

For the total ground state energy as a function of the radius of the shell we 
get the curves shown below (Fig. 3.3-3.5) for different values of the strength of 
the potential a. 




Figure 3.3: The renormalized vacuum energy 
values of the potential. 



ren 
9 



[R) for positive and negative 
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Figure 3.4: The renormalized vacuum energy E"^"'{R) for a=0.3 



3.7 Discussion 



We have obtained a representation of the renormahzed ground state energy of 
a scalar massive field in the background of a semi-transparent shell containing 
convergent integrals of simple functions. This expression is given by the sum 
of (|3.49|) and ( |3.51| ) and it depends only on the two parameters of the classical 
system, namely the radius and the strength of the potential of the spherical shell. 
The plots of E^^"' as a function of the radius show that for repulsive potentials 
the renormalized ground state energy is positive only in some limited intervals of 
the radius axis and only when a is smaller than 1. For a strength of the potential 
larger than 1 the energy is always negative. This is the most striking conclusion 
of this chapter. For very large values of a the shell should be no more transparent 
and the problem should formally become a Dirichlet boundary condition problem. 
One could check this in the equation ( p.32|) for the Jost function: here inserting 
a large a the addend 1 becomes negligible, then one would just have the product 
of the two modified Bessel I and K functions; such a Jost function is exactly the 
one for a perfectly reflecting spherical shell (Dirichlet boundary conditions). In 
that case the ground state energy is simply the sum of the energies inside and 
outside the shell. Then formally we should have: 



lim 

a— >+oo 



GSE 



semi— trans. 



GSE' 



(3.52) 



Now it is shown in that the "mirror" configuration, in the massive field 
case, has always positive ground state energy for repulsive potentials. This is in 
contradiction with our plots which show an opposite sign. Furthermore, the A2 
coefficient in paper is zero. In our work A2 remains a non zero coefficient also 
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Figure 3.5: The renormalized vacuum energy E^^"'[R) for a=10 



in the hmit a — > +00 demonstrating that the transition hypothesized in ( |3.52| ) 
is singular. For flat parallel semi transparent boundaries with delta function 
potential in the vacuum of a scalar massive field, the transition is actually fulfilled 
as shown in paper [|20|. In the configuration analyzed in this paper the limit 



( |3.52D works only for the regularized GSE, as mentioned above, but after the 
renormalization the limit is no more valid. This means that the subtraction of 
the divergent part of the energy and the limit a —>■ +00 are two non commutative 
operations. We remark again that the results of this work cannot be directly 
applied to the case of a massless field, since the initial normalization condition 
would fail and the vacuum energy would not be univocally defined. 
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Chapter 4 



Vacuum energy for a 
semi-transparent cylindrical shell 



Calculations with cylindrical geometries began historically with |^9[, after the 
suggestion that a cylinder, as a kind of intermediate shape between the parallel 
plates and the sphere, could possess a zero Casimir stress [^. However paper 



showed that a perfectly conducting cylinder in the electromagnetic vacuum 
has a negative Casimir energy. A number of papers on the dielectric cylinder 



3^ , have recently showed the interesting result of a vanishing vacuum energy in 
the dilute case. We find of interest to carry on this research with the analysis 
of massive fields in the background of non-ideal cylindrical boundaries. In this 
chapter we consider a hollow cylinder with radius R, having a delta function S{r — 
R) as a potential profile. This model can be considered as a "scalar" version of the 
dielectric background. A non singular potential would be, of course, more realistic 
under the physical point of view, but the calculations would be considerably 
more complicated. We hope this work can contribute in the understanding of 
the structure and meaning of vacuum. The cylindrical shell could, besides, find 
physical applications in the calculation of the quantum corrections to vortices in 
QCD or in the electroweak theory. Another interesting perspective is in the recent 



discovery of the so called nanotubes [Q, which are large carbon molecules 



generated in laboratory offering the intriguing possibility of measuring quantum 
effects on small cylindrical objects. 

The setup of the problem is analogous to that of the spherical shell. We 
will first approach the quantum field theoretical problem and discuss the renor- 
malization. Then, the Jost function will be calculated. In the remaining part 
of the chapter, the numerical computation will be performed and graphically 
represented. 
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4.1 Vacuum energy in terms of the Jost function 
and renormalization 



Let us work with a real massive scalar field (f){x) with mass m and let us quantize 
it in the background of a cylindrical potential. The field equation in cylindrical 
coordinates r, (p, z, after separation of the variables, reads 



pl-m''-pl--- V{r) + -dr + dl ) 0Kpo,Pz, r) = , (4.1) 



1 

where is the momentum four vector being pz its component along the longi- 
tudinal axis of the cylinder and / is the angular momentum quantum number. 
V{r-) is the background potential given by 

V{r) = ^S{r - R) , (4.2) 

it represents an infinitely thin cylindrical shell whose profile is a delta function. 
The shell has a circular section of radius R and it extends from — oo to +00 in the 
z direction, a is the dimensionless parameter giving the strength of the potential. 

Considering eq.(4.1) as a scattering problem, we choose the "regular solution" 
which is given by 

<Pir) = Ji{kr)Q{R-r) + ^{fi{k)Hl'\kr) + f;ik)Hl'\kr))e{r-R), (4.3) 

where k — yjp\ — w? — pi, Jm{kr) is a Bessel function of the first kind, Hi^\kr) 
and H\^\kr) are the Hankel functions of the first and of the second kind and 

the coefficients fi{k) and f*{k) are the Jost function and its complex conjugate 
respectively. 0(i? — r) is a theta function. The field is therefore free in the regions 
< r < i? and R < r < 00. At r = R the field is continuous and we have the 
following matching conditions 

0'(i? + O)-0'(^-O)= |0(i?) 

(j){R + 0) = (PiR - 0) , ^ ^ ^ 

where the prime indicates derivative with respect to r. 

The quantum field is in its vacuum state, its energy is given by half of the sum 
over all possible eigenvalues ur^n) of the Hamilton operator related with the wave 
equation (4.1). We define a regularized vacuum energy 

(n) 

where s is the regularization parameter which we will put to zero after the renor- 
malization, /I is the mass parameter introduced for dimensional reasons and (n) 
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includes all possible quantum numbers. We will calculate the energy density per 
unit length of the cylinder given by 

°^ (n) 

where e(„) are the eigenvalues of the operator contained in (4.1) without pz and 
with k = -s/pq— m?- We carry out the integration over pz and we arrive at 

Following a procedure described in the last chapter^ we transform the sum in 
(4.7) into a contour integral and, dropping the Minkowski space contribution, we 
arrive at 

£ = —Cs / dkie + mY~'dk\nfi{tk) , (4.8) 

l=—oo 

where Cg = (1 + s(— 1 + 2 ln(2/i)))/(27r) and fi{ik) is the Jost function defined in 
(4.3) on the imaginary axis. It contains all the information about the background 
potential under examination. We will find fi{ik) explicitly in the following sec- 
tion. The energy defined in (4.8) is renormalized by direct subtraction of its 
divergent part 

£ren £ £div (^'9) 

with the normalization condition demanding that the vacuum fluctuations vanish 
for a field of infinite mass 

lim £ren = . (4.10) 

m— >oo 

This normalization condition eliminates the arbitrariness of the mass parameter 



/i (the interested reader can find more details in |^). The isolation of £div is 



done with the use the heat-kernel expansion of the ground state energy 

r(s + j-2) 



£ = y V^ m^-^(-+^-)A, , J = 0, ^, 1, ... , (4.11) 

^ 327r2 T s + l 1 ^ J '2' ' ' ^ ' 

j 

where the Aj are the heat kernel coefficients related to the background. Then it 
is possible to define 

m? f 1 , 4u^ \ ^ "m ^ 

327r^ \s / IbiT'^i^ 

—J- + \n^-2\A2, (4.12) 



^ The procedure is explained with details in 
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which includes all the pole terms and all the terms proportional to non-negative 
powers of the mass. 

In order to perform the analytical continuation s — * we split the renormal- 
ized vacuum energy (4.9) into two parts 

£ren = £f + ^as, (4-13) 

with 

^/ = -i^^ E y dk[e-m'Y'^ — [\nm)-\nfrm] (4.14) 

and 

£as = —^C^Y, j dk[k^~m^Y-^— In frm-Sd.,. (4.15) 



/=— oo 



Here fj^^{ik) is a portion of the uniform asymptotic expansion of the Jost function 
which must include as many orders in / as it is necessary to have 

In fiitk)- In =0{l-') (4.16) 

in the limit I —>■ oo, k oo with fixed l/k, which is sufficient to let the sum 
and the integral in (4.14) converge in the limit s = 0. The splitting proposed in 
(4.13) leaves the quantity Sren unchanged, while it permits the substitution s = 
in the finite part Sf. 



4.2 The Jost function and its asymptotics 



We insert solution (4.3) into (4.4) and we find 



Ji{kR) = i fi{k)HP{kR) + ft{k)H\'\kR) 



(1), 



fl{k)Hf\kR) + ft{k)H\'\kR)^) \r=R = fJlikR) + {drJlikr)) \r=R . 

(4.17) 

The system can be solved for fi{k), with the help of the Wronskian determinant 
of the Hankel functions [^|, the result is 



r(l)/ 



m 



-UnkR {^Ji+i{kR)Hl^\kR) - Ji{kR)H}l\{kR) 



a 



-—J,(kR)HrikR) 



rW(UG\\ (4-18) 
The corresponding Jost function on the imaginary axis can be written in terms 



of the modified Bessel I and K functions, again with the help of ||5^ 

Mik) = l + a IiikR)Ki{kR). (4.19) 
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From the Jost function (4.19) one arrives at the uniform asymptotic expansions 
fi'^^iik) for positive / and fi^~{ik) for negative / and at the asymptotic expansion 
fo^{ik) for I = 0, with the help of the asymptotics of the Bessel I and K functions 
for large indices and large arguments available on |^ . Since the asymptotic Jost 
function consists of three different contributions, the sum over / appearing in 

(4.14) and (4.15) must also be distinguished in three contributions: a sum over 
negative I, a sum over positive / and a contribution coming from / = 0. The first 
two contributions can be summed up analytically in the following way 

— 1 oo oo 

...ln/--(^A;)... + 5^...1n/-+(^A;)... = ...(ln/-+(^fc) + ln/--(^fc))... , 

l=-oo 1=1 1=1 

(4.20) 

where the dots represent, for simplicity, the rest of the functions in (4.14) and 

(4.15) . Then, eq.(4.14) and (4.15) are rewritten in the form 

1 CO „oo o 

£as = -lCsJ2 dk[k'-m']'-^—\nfr^itk) 
1=1 

—Cj dk[k''-m'Y-^—\nf^'{lk)-£a^v (4.21) 
4 dk 

and 

= -^Y.]^ dk[e - m']^ {2ln Mtk)- In fr^{^k)) 

-— y dk[e-m']— {In fom-\nfS'{zk)) , (4.22) 

where In Z"**^ (ifc) = In f^^'^'^ (ik) + ln/";~(iA;) and we have used the property 
fi{ik) = f.i{ik) of eq.(4.19). 

Taking the logarithm of the uniform asymptotics of the modified Bessel func- 
tions and re-expanding in negative powers of the variable I (see for details 
on this procedure, see also appendix A) we find 

ln/<"±(ifc) = X^J^X,,^^, (4.24) 

n=l j 

where t = (1 + {kR)/lY)2 and the non-vanishing coefficients are 

Xi,i = a , X2,2 = -0^/4 , 

^3,3 = «/8 + aVl2 , X3,5 = -3a/4 , (4.25) 
Xgj = 5a/8 . 
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In this definition we have included 3 orders in I which are sufficient to satisfy 
condition (4.16). Substituting (4.23) and (4.24) in (4.21) and (4.22) we find 



I — 1 \^ — 1 j J 



and 

£f = 



__L £ - (ln/„(*) - ^ - ^) . (4.27) 

We call the first addend in (4.27) Sji and the second Efo, that is Sf = Sfi + SfQ. 

4.3 Asymptotic part of the energy 

We go forward with an analytical simplification of (4.26). We call the second 
addend in (4.26) Saso, it can be immediately calculated, giving 

A simplification of the first addend in (4.26) can be achieved with the Abel-Plana 
formula 



in the present case the function F(/) is 



/•oo / 3 

F{1)=\ dk{k^ + m^)'-'dAY.Y.^- 

■^'^ \n=l 1 



3 



(4.30) 



Then, the first addend in (4.26) turns out to be the sum of three contributions. 
The first contribution can be calculated by means of formula (B.l) displayed in 
the appendix. We find 
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The integration of the second term in the Abel-Plana formula can be performed 
with formula (B.2), which yields immediately 

-2 /I , /v^ 



as2 



am a 
+ 



SnR 647ri?2 V s 



+ In 



a 



967imR^ 

The third contribution, after application of formula (B.2), reads 

1 ''^ 

Sas3 = --C,^X„,,(-m^2T(2-s))Z,,,(mi?) 



(4.32) 



(4.33) 



where 



Zn,j{mR) 



dl P-"-2sin[f(j-n)] 



1-e 



+ 



mR 



dl p-"-2sin[7r(l-5-f)] 



.+ §-1 



(4.34) 



Here, the integration over / cannot be performed analytically. However the ex- 
pression can be remarkably simplified by integrating several times by parts. We 
find finally 

2 

rv r' 



hi{mR) — — — - h2{mR) 



27ri?2 



32i?2 



+ 
+ 



a a'' 
+ 



167ri?2 247ri?2 
a 



h^{mR) — 



a 



/i4(mi?) 



hr,{mR) , 



(4.35) 



487ri?2 

where the functions hi{x) are convergent integrals analogous to those found in 
the last chapter 

""^ dl 



hi{x 
h2{x 

h^{x 
h^{x 



1-e 
dl 

dl 



2ttI 



X^ 



1 - e2-' / 



^7 d'--') 



2 _ J.2 



dl 



dl 



I 



1 - 62^^' 
/3 



2 _ rj.2 



l-e2^W // I 



2 _ rf.2 



(4.36) 
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In appendix B the reader can find more details about the derivation of (4.31), 
(4.32), (4.35). The contributions Sasi and £as2 contain all the pole terms (all the 
divergences of the vacuum energy) plus the terms proportional to non-negative 
powers of the mass (which do not satisfy the normalization condition). All these 
terms are subtracted and are used to calculate the heat-kernel coefficients by 
means of definition (4.12) for Sdiv Below we give the heat kernel coefficients which 
we calculated up to the coefficient (adding four more orders in In f^"^^ (ik)) , in 
the hope that they will be of use for future investigations on the same background 



2^3/2 



a TV' 



Ao = , = 

Ai = -27ia , A3/2 - 

4 _ TTQ^ A - (3a^+4Q:4)^3/2 

A _ (4a^+7a^)TT a _ (81o2+l20a''+128a'')7r3/2 

— mnp4 , ^7 



^3 ' 7/2 24576iJ5 

(640=^+5205+390-7)^ 



A , — ^ 

4 16380/?6 

One can note how all the integer heat-kernel coefficients depend on odd powers of 
the coupling constant, while the half- integer coefficients depend on even powers 
of a. The same feature is present in the heat kernel coefficients of a (5-potential 
spherical shell (section 3.4 of the preceding chapter). We note also that (4.32) 
is in agreement with the coefficients that one would obtain from lemma 2 of 
paper [Q, where the heat kernel expansion for semi-transparent boundaries in 
(i-dimensions is examined. 

We perform the subtraction of the divergent portion and we obtain the final 
result 



a. a 
^as = :t hiiniR) — hoimR) 

+ I — H I himR) 

OL OL 01 

-8^^^^"^^) +48^^^^"^^) - 64^ 

4.4 Finite part of the energy and numerical re- 
sults 

The quantity (4.27) cannot be analytically simplified. We integrate £fi and £fQ 
by parts and make the substitution k —>■ k/R to get an explicit dependence on R 



" 4;;^££^^^ ^ (^21n/K«/^)|fc^fc/fl- ) 1^^^^^^ 
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For small values of R, Sfi behaves like R~^, while SfQ has a logarithmic behaviour 

a^ln(mi?) /, 

We note that in the limit a 0, the logarithm of the Jost function fi{ik) can be 
expanded in powers of a, giving (from eq.(4.19)) 

\n{l + aIi(kR)Ki{kR)) aIi{kR)Ki{kR) + ^a^lf{kR)Kf{kR) + 0{a^), 

then, the summation over / of the leading term of this expansion could be an- 
alytically performed following a method recently proposed in [^. This would 
give us a first order approximation of Sfi. However we will restrict us here to a 
numerical calculation of the sum in (4.39), which can be performed for small as 
well as for large values of a. 

To find the asymptotic behaviour of £as we rewrite eq.(4.38) in the following 
form 

Sas = [awi{mR) + a^W2{mR) + a^W3(mi?)] , (4.41) 

where the functions wi(x), W2{x) and w^{x) are given by 

/ 111 

wi{x) = [ -hi{x) + -h3{x) - -h4{x) + —h5{x) 



4 -V / 24 ' ' 32x 

W2{x) = [-Y^h2{x] 

= (^ft3(.) - ^) . (4.42) 

As mentioned above, the functions hi{x) are quickly converging integrals. The 
behaviour of £as is governed by the functions ^1,2,3(2;) and by the value of a. For 
i? — > we find 

wi{x) ~ 0.0000868 + 0{x) , 

W2ix) ~ —lnx + 0.115 + 0{x) , 
16 

wsix) ~ ^Inx + 0.00479 + 0{x) . (4.43) 

Thus, for small values of the radius, Sas is proportional to Ini?. Summing the 
contribution coming from (4.40) and (4.43), we find that the leading term of the 
renormalized energy, for i? — 0, is 

Sren ■ 4.44 
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This result is in agreement with the prediction (section II of paper p9| ) 



liraSren ~ In (mi?) . (4.45) 

In the limit i? — oo the behaviour of the renormalized energy is determined by 
the first non- vanishing heat kernel coefficient after A2. From (4.37) and eq.(4.11) 
we arrive at 

lim Eren ^ ^ + O — . (4.46) 

ij^oo 2048mi?3 1536mi?3 ^ yR'^ ) ^ ' 

We have numerically evaluated the quantities £as^ £fh ^fo and £ren as func- 
tions of R, fixing the values of the mass to 1. It turned out to be necessary 
to sum 20 terms in the variable / and to integrate up to 1000 in the variable k 
to obtain "stable" numerical values for the energy. Below we give the plots of 
the various contributions to the vacuum energy and the complete renormalized 
energy for different values of the potential strength. For a < we found the 
renormalized vacuum energy to posses a small imaginary part becoming larger 
when a approaches —oo. Particle creation accounts for this contribution. It 
starts when the attractive potential of the shell becomes over-critical, that is 
when e < —m < 0, where e is the energy of the bound state; in this case the 
effective action of the system acquires an imaginary part0, however a detailed 
discussion of this aspect of the theory is beyond the purpose of the present chap- 
ter. It should only be said that in the plots traced for negative values of a (Fig. 
4.4 and 4.5) the energy is to be intended as real part of. 



4.5 Discussion 

In this chapter we calculated the vacuum energy of a scalar field in the background 
of a cylindrical semi-transparent shell. The formulas for the energy density per 
unit length of the shell are given by equations (4.38) and (4.39). The heat kernel 
coefficients (4.37) for a cylindrical potential containing a delta function, are also 
a relevant part of the results. A discussion of the sign of the vacuum energy is 
possible by means of the data given in the numerical section of this chapter. The 
energy is found to be negative in the background of repulsive potentials (a > 0) 
for every finite value of the radius (Fig. 4.3). This conclusion sets a closeness 
between the model studied here and that of a conducting cylinder |^ and of a 



dielectric cylinder |5^ in the electromagnetic vacuum. The latter models possess 
also a negative Casimir energy. There is also a resemblance with the ^-potential 



^For the theoretical foundations of this phenomenon see IQ. In the context of vacuum 
energy the appearance of an imaginary part in the renormaUzed energy was observed for instance 
in0. 
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Figure 4.1: Repulsive potential. The contributions to the renormalized vacuum 
energy multiplied by i?^ • o;"^, for a = 2.1 . 
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Figure 4.3: Repulsive potential. The complete renormalized vacuum energy 
8ren{R) multiplied by E? ■ o;"^, for a = 0.3, a = 2.1 and a = 10. 
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spherical shell (see last chapter) investigated in paper where a negative 
energy was observed for large repulsive potentials. 

In the background of attractive potentials {a < 0) the energy is negative on 
almost all the R-axis. As we see in Fig. 4.5, the energy becomes positive for very 
small values of the radius. This happens when the logarithmic term in (4.44) 
compensates for the term proportional to — in (4.46) which is dominant for 
small values of a in the region R> 1. The thin region of positive asymptotically 
increasing energy appearing in Fig. 4.5 is reasonably out of physical applicability. 
For both repulsive and attractive potentials the renormalized vacuum energy of 
the semi-transparent cylinder goes logarithmically to ±oo in the limit R ^ 0. 
This logarithmic behaviour was also found in the semi-transparent spherical shell, 
however it was not found in the 5-potential flux tube |Q (see chapters 5 and 6), 
which has many feature in common with the background investigated here. This 
can be explained by means of the heat-kernel coefficient A2, which is a non zero 
coefficient here as well as in ||6^. It is also interesting to note how in expression 
(4.37) for the coefficient A2 the contributions proportional to a and have 
cancelled and only a term proportional to is present. The cancellation of the 
lower powers of the potential strength was also observed in [^] and in paper 
where a theta function profile in a dielectric spherical shell is examined. 
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It confirms the observation that singular profiles show weaker divergences than 
smoother, less singular profiles. 



R 



Figure 4.5: Attractive potential. The complete renormalized vacuum energy 



8ren{R) multiplied by R^ 



■ a 



for a = —0.3. The peak beneath the lowest point 



of the curve is due to the presence of a small imaginary part in the renormalized 
energy causing the function to be no more analytical. The energy starts to be 
complex at -R ~ 0.04. 
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Chapter 5 



Vacuum energy in the presence of 
a magnetic flux tube: Scalar case 



The formalism applied in the last two chapters in the context of curved mate- 
rial boundaries, can be extended to problems with special background fields, like 
classical magnetic strings. The calculations will follow as usual, provided the Jost 



function is found. This method was used in paper to analyze the vacuum of a 
spinor field in the presence of a homogeneous magnetic flux tube of finite radius. 
The vacuum energy was found to be negative and it did not show a minimum 
for any finite value of the radius. A natural question is if inhomogeneous mag- 
netic fields can minimize the energy and render the string stable. The question 
was already raised in [|^. The present paper extends the investigation begun in 



5^ to an inhomogeneous magnetic string with delta function profile. The delta 
function, represents here a simple example of inhomogeneity, which could give 
an insight into the problem of vacuum fluctuations in magnetic backgrounds. It 
is a kind of "semi-transparent magnetic boundary" . The quantum mechanics of 
spinor fields in the presence of magnetic fluxes has been elaborated in early works 



33],p5[ while more recent investigation in this direction has been motivated by 



the interest in the Aharonov-Bohm effect ||36|, Singular inhomogeneous mag- 
netic fields were examined in for the calculation of the fermion determinant 
and in for the investigation of the bound states of an electron, however the 



ground state energy was not calculated in those works. In this chapter the ground 
state energy will be calculated for a scalar field, while the spinor field will be con- 
sidered in the next chapter. In the first part of this chapter we calculate the Jost 
functions and the heat-kernel coefficients. The energy will be renormalized as for 
the preceding chapters, imposing the vanishing of the vacuum fiuctuations for 
fields of infinite mass. In the second part of the chapter we will work numerically 
on the renormalized energy finding its asymptotic behaviour for small and for 
large values of the radius of the string. 
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5.1 Scalar electrodynamics 



We quantize a scalar field $ in the presence of a classical magnetic field whose 
form is that of a cylindrical shell with delta function profile. The section of the 
string is a circle with radius R. The magnetic field is given by 

B{r) - l^Kr- R) (5.1) 

where is the magnetic flux, r = \/ x"^ + y"^ and z is the axis along which the 
cyhndrical shell extends to infinity. The quantum field obeys the Klein-Gordon 
equation for the scalar electrodynamics: 

(D^ + m^)$(a;) = (5.2) 
where mg is the mass of the field, and 

D2 = 52 _ v2 - 2ieA^dt - 2ieAV - e^A^ ^ + e^A^ , (5.3) 

here e is the electron charge, is the vector potential of the electromagnetic 
field and the convention g^^ =diag(l, —1, —1, —1) is used as well as the gauge 
d'^A^j^ — 0. We want to find a solution of eq.(5.2) in cyhndrical coordinates z,r, 
</?, then we take the relevant operators in (5.3) in the following form 

•f^ ( sina?^ . ^ cos(z3^ ^\ ,^ 
V — > I — cos</7ar a<^, sm (^a^H a<^, j (5.4) 

V^- \rd,^^^dl^dl. (5.5) 

The potential four vector A^ associated with the magnetic field (5.1) contains a 
theta-function: 

^^±_Qir^ ^" = 0; (6.6) 

27r r 

where = (— sin 0, cos0, 0). Therefore operator (5.3) becomes 

Here, and in the rest of this chapter and of chapter 6, P will represent the strength 
of the background potential 

/J = I* . (5.8) 

With the ansatz of the separation of the variables the scalar field is transformed 
into 

$(x) exp{ipot-ip^z + imip) ^rn{po,Pz,r) , (5.9) 
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where is the momentum four vector and m is the orbital momentum quantum 
number. Combining (5.2), (5.7) and (5.9) the new field equation reads 

Pl-ml-pl- ^"^' " + ^dr + d^'^ $^(po,P.,r) =0. (5.10) 

or 

t^- - ''"'^^'r^»' + ia. + 9;) *„(/l. r) = . (5,11) 

where k = \/pq — ml — p^. The solutions of this equation are Bessel and Neu- 
mann functions. The kind of function and their coefficients are to be determined 
by means of physical considerations. We take here the regular solution which in 
general scattering theory has the known asymptotic behaviour 

$ '■~°J^(fcr), (5.12) 



$ ™ ^iumi^Ukr) + nmi^i.ikr)), (5.13) 

where Jm{kr) is a Bessel function of the first kind, Hl^^_ij{kr) and H'!^_p{kr) are 
Hankel functions of the first and of the second kind and the coefficients fm{k) 
and f^{k) are a Jost function and its complex conjugate respectively. In the case 
of the delta function profile the regular solution reads 



$(r) = J^(fcr)e(i?-r) + i(/„(A:)ifi%(A:r)+/;(A;)iJ«^(A:r))e(r-i?) , (5.14) 

then, we can define a field in the region r < R inside the cylinder 

^l{k,r) = J^kr) (5.15) 

which is independent of the strength (3 of the potential, and a field in the 
region outside the cylinder r > R 

^Zik,r) = l{Uk)H^lp{kr) + nk)H^J,lp{kr^^ (5.16) 

which describes incoming and outgoing cylindrical waves. The conditions for the 
field at r = will be discussed later. 



5.2 Ground state energy in terms of the Jost 
function and normalization condition 

A regularized vacuum energy can be defined as 
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where the e(^n,a) are the eigenvalues of the Hamiltonian operator associated with 
(5.2), a = ±1 being the index for the particle-anti-particle degree of freedom, 
while n includes all oder quantum numbers, s is the regularization parameter to 
be put to zero after the renormalization and is the mass parameter necessary 
to maintain the correct dimensions of the energy. The string is invariant under 
translations along the z axis, therefore the energy density per unit length of the 
string is 

= l^'^r ^ T.ipi+>^uy''~' ' (5-18) 

°° (n,a) 

where the A(„) are the eigenvalues of the operator defined in (12) with k — 
•\/Pq — ml . We perform the integration over pz in (5.18), getting 

{n,a) 

The next step is to transform the summation in (5.19) into a contour integral. We 
enclose temporarily the system in a large cylindrical quantization box imposing 
Dirichlet boundary conditions for the field at this boundary like we did in section 
3.2. The final result will be independent of this choice, since the quantization box 
must be expanded to infinity. Now we can express the eigenvalues A(„) by means 
of the zeros of solution (5.16), which is valid in a region far from the origin. The 
derivative of the logarithm of the sohition will then have poles at A; = A(„). We 
can express (5.19) through an integral whose contour encloses these poles which 
lie on the real axis k. The deformation of the contour on the imaginary axis and 
the dropping of the Minkowski space contribution allows to reach the final form. 

= -oC. E / dk{k^ + miy-^dk\nU{ik) , (5.20) 

m=-oo "^"^e 

where fmiik) is the Jost function with imaginary argument and Cg = (l + s(— 1-|- 
2 ln(2/i)))/(27r) is, as in chapter 4, a simple function of the regularization pa- 
rameter. The renormalization of (5.20) is carried out by direct subtraction of its 
divergent part 

S^en = - ■ (5.21) 

The isolation of the divergent part from the total energy will be performed via 
heat-kernel expansion as we will see in a moment. The subtracted part should be 
added in the classical part of the energy resulting in a renormalization of the clas- 
sical parameters of the string, however we do not treat the classical energy of the 
system here but only the vacuum contribution. For the analytical continuation 
s — > we split £^g„ it into a "finite" and an "asymptotic" part. 

= + (5.22) 
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with 



and 



1 /-oo 

^r--2^s E / dk[e-ml]'-^ — [In Uik)- In f^^ik)] (5.23) 

m=-oo "^"^e 

^'at^-^Cs X, ^^[^' - ^e]'-^^ In /-(i^) (5.24) 



m=— oo 



where /^f is a portion of the uniform asymptotic expansion of the Jost function. 
The number of orders to be included in must be sufficient to let the function 

ln/^(iA;)-lnC(iA;) (5.25) 

fall as (or /s"^) for k and m equally large, in this case the integral and the 
summation in (5.23) converge for s — > 0. To this purpose the usual three orders 
in the asymptotics are enough. The splitting proposed in (5.22) immediately 
permits the analytical continuation s = in S'p, furthermore it allows a very 
quick subtraction of the pole terms in the asymptotic part (5.24). S^g is a finite 
quantity. For the definition of fjf^ we write the sum (5.19) as 



2s 



2r(s- 1/2) 



/•oo 

/ dt t'-^/^K{t) (5.26) 
Jo 



where K(t) is the heat kernel related to the Hamilton operator , which can be 
expanded for small t 

_(^2 OO _ 

m = - 7|^E^.■^^ j = o,-,i... (5.27) 

The Aj are the heat-kernel coefficients related with the Hamiltonian operator. 
By means of (5.26) and (5.27) we can expand the ground state energy in powers 
of the mass and get 

CSC _ V ^('^+^'-2) 4_2(s+,-) , . 

j ^ ' 

in which the divergent contribution is isolated 

A /I ^,,2 



CSC 



ml 



1 1 ^Ai^ 1\ ml . 



647r2 
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The poles are all contained in the three terms corresponding to the heat kernel 
coefficients ^0,^1,^2, however we included in E^^^^ two more terms in order to 
satisfy the normalization condition 



lim Ell^ = . (5.30) 
Again, the normalization condition must be changed if one desires to investigate 



the limit m — 0. In paper |3T[ the interested reader can find further details and 
comments about this limit. 



5.3 The Jost function and its asymptotics 

We calculate now the Jost function fm{k) related with the delta function potential 
of the magnetic string. We have defined a field in the outer region r > R and a 
field in the inner region r < R, we now need some matching conditions on the 
boundary r = R. As the delta function is a continuous function, we require that 
the field is continuous on the boundary. From this condition and from the field 
equation (5.11) it follows directly that the ffist derivative of the field must be 
continuous on the boundary: 



and with the use of solution (5.15) and (5.16) 

drJm{kr)\r=R = ^ (/n^ ( A:)i/i% ( A:r ) + /;;(fc)i/i% ( A;r ) ) | 
Solving for fmik) one finds 



(5.31) 



(5.32) 



^ mJm{kr)\r=RH^^l^{kR) - drH^^\p{kr%=RMkR)) 
H;^UkR)drHi^l,{kr)\r=R - H^Zp{kR)dMl,{kT)\r=n ' 



which with the use of the Wronskian determinant for the Hankel functions p6 
becomes 



- Jm^i{kR)H^^l,{kR) + ^Jm{kR)H^^\p{ 



(5.34) 

The corresponding Jost function on the imaginary axis is written in terms of 
modified Bessel I and Bessel K functions 

]+iP(3ImK^.p . (5.35) 
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where the arguments (kR) of the Bessel functions are omitted for simphcity. Eq. 
(5.35) can be written in a more compact form 

U{ik) = i^kR[l'^Km-p + ImK'^.p\ . (5.36) 

where the prime denotes the derivative with respect to the argument. This ex- 
pression holds for positive and for negative values of m. On the contrary the 
uniform asymptotic expansion, which we need for formulas (5.23) and (5.24), is 
a different function for positive or for negative m. To find it we write the Bessel 
I and K functions of eq.(5.35) in the form 

Im+a{,{rn + a)z'), Km+a{{m + a)z'), z' = kR/{m + a), (5.37) 

where a can be or 1 for the Bessel I function and —jS or —f3 + 1 for the Bessel 
K function. Their expansions for large positive orders are well known [Q, for 
instance Km+aii^ + a)z') is expanded as 



y 2(m + a) (1 + z'^)^'* y ^—^ (m + 



(5.38) 

where r] = Vl + ^'"^ '^'^'^{ ij^Jij^^n ) ^^"^ t^e Debye polynomials in the 

variable x. However we are interested in an expansion in powers of the variable 
m alone: an expansion of the form 

K^U{^ + (^)z) ~ , (5.39) 

n 

where the X„ are some coefficients depending on k, R and /3. Therefore we make 
the substitution z' ^ z {:;^^), with z = KR/m in the argument of the Bessel 
function and we rewrite its expansion as 




-(m + a)T? 



with the obvious change for rj' . Then, re-expanding in powers 1/m" we get 



3 



Km+a{kR) ~ W — exp{> m-"^^(n,a,t)} , (5.41) 
V 2m ^-^ 

n=— 1 

where t = (1 -|- 2;^)"^/^ and the functions 5'K(n, a, t) are given explicitly in the 
Appendix. The corresponding expansion of the Bessel I function in negative 
powers of m will be 

1 ^ 

lm^a{kR) ~ -^==exp{ J] m-"5/(n,a,t)} , (5.42) 
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where the functions Si{n,a,t) are given in the appendix. Inserting these expan- 
sions in (5.35) one finds an asymptotic Jost function vahd for positive m, we 
name it f^~^{ik). To find the asymptotics for negative m we must take the Jost 
function in its form 

fm{ik) = i^kR [I_rnK^rn+(3-l + I-m-lK-m+p] + l'^ P I-mK -m+l3 , (5.43) 

which is identical to (5.35) because of the property of the modified Bessel func- 
tions 

Iiix) = I^i{x) , K^{x) = I^^{x) , (5.44) 

where / is any natural number and p any real number. Then in (5.43) a large 
positive index always correspond to a large negative m and inserting (5.41) and 
(5.42) in (5.43) we find an asymptotic Jost function valid for negative m, we call 
it /^-(^fc). 

For m = 0, the asymptotic Jost function can be obtained from (5.35), using 



the expansions |56 



we call this contribution /q'*(?A;). The finite and the asymptotic part of the energy 
defined in (5. 23), (5. 24) are also split into three contributions: one for positive m, 
one for negative m and one for m = 0. The positive and negative contributions 
can be summed up in a single term, but the contribution coming from m = 
must be calculated separately and summed just numerically at the end, in fact 
we have 



/ oo „c 



= "2^^ II./ dk{e + ml)d,{\nf^{zk)-\nf:^^^{zk)) 



oo 

2 I _2 



and 



where 



+ / dk {e + mi)dki\nfoitk)-\nfS'itk)) , ) (5.47) 

1 / °° poo 

^at = -1^CJJ2 dk{k' + mlY-^^,\nf:^^{^k) 

+ rdkie + miy-'d,lnfS^{tk) - S^}j , (5.48) 

J me J 

In = In frn{ik) + In f-miik) 
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and 

ln/-±(zfc) =ln/-+(zfc)+ln/-^-(zfc). 

Taking the logarithm of f!^^{ik) and /"^(i/c) and re-expanding in powers of m 
we find the functions needed in the integrand of (5.47) (5.48) up to the desired 
order. We are interested in In f^^{ik) and \nfQ^{ik) up to the third order, they 
read 

ln/„" = ^ , (5.49) 

n t 

where t = 1/(1 + (kR/mY)^ and the nonzero coefficients are 
^1,1 = , ^2,4 = /3V4 , 

X3,3 = /3V24 - PVl2 , X3,5 = + /574 , (5.51) 

X3,7 = /3Vl6 . 

As we found also in the preceding chapters, three orders are sufficient for the 
convergence of S^^. More orders could have been included in definition (5.49) 
and (5.50) to let the sum and the integral in (5.47) converge more rapidlyQ. 



5.4 The asymptotic part of the energy and the 
heat kernel coefficients 

Having found the Jost function related to the cylindrical delta potential an impor- 
tant part of the calculation is done. We proceed with the analytical simplification 
of S^g. The second term in (5.48) which we name S^'^q, can be quickly calculated 
inserting in it the result (5.49). We find 

Co - -'^ ■ (5.52) 

The first term in (5.48), which we name here ^^^^(m) contains the sum over m 
which can be carried out with the Abel-Plana formula 



5^F(m) 



m=l 



, , 1 ^/ N dm Fiim) — Fi—im) 
dmF{m) - -F(0)+ / ^ _ ^ \. ^ ■ (5.53) 



^ We would like to stress that with the introduction of asymptotic expansions in our calcu- 
lation we do not approximate the vacuum energy. The total energy as defined in (5.21) remains 
an exact quantity. The uniform asymptotics of the Jost function is just a mathematical tool 
which permits the analytical continuation s — > 0. 
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In our case is 



F{m)= dk{e + miy-'dklnf^'^{ik) . (5.54) 

Jnie 

Thus is split into three addenda: 

1 /"OO /"OO 

^Tsi=-7^Cs dm dk{k' + ml)'-'dklnf^^{tk) (5.55) 

^ -/o Jme 

^ts2 = ^^^.i^(0) , (5.56) 
CSC _ V r dm F{im)-F{-im) 

The contributions and £^^2 can be calculated with the formulas given in 
appendix B. The result is 



e 



The divergences are all contained in the first term of ^^^si- "^^^ fi^^^t term of £^^2 
will cancel with f^^^Q and it remains only one term containing a positive power of 
the mass. This term is the contribution to Ef^ corresponding to the heat kernel 
coefficient ^3/2 (see eq.(5.29)) and therefore it will be subtracted as well as the 
pole term. 

The last contribution to is £^^^3, whose calculation demands a little more 
work. Using the formula displayed in appendix B to calculate the integral over 
/c, we we find 

^asZ = (-l)^^^.E^"-^- [-^e-'^r(2 - S)\ KA^eR) , (5.60) 

where the functions Anj{meR) are given by 
A„,(x) ^ ^>±^ 



fX dm " 2sm[f (i-n)] 

Jo l-e^^rn T 



1-2 



I roo dm m-'-" 2 sm[7r(l-s-n/2)] 



(5.61) 



We have calculated these functions by partial integration for n < 3 and j < 7 
and we show them explicitly in the appendix B. By means of those functions and 
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of the coefficients (5.51) we find 




' 1 



m 



rUeR) 



2 



16R 1 - e2'^"^=« ' 



(5.62) 



where the prime in the integrand denotes derivative with respect to m. The last 
term in (5.62) goes to zero for rrie oo and therefore must be included in the 
renormalized energy, however, given its exponential behaviour, it does not con- 
tribute to the heat kernel coefficients. The heat kernel coefficients, which we have 
calculated up to the coefficient ^47/2 (including four more orders in In (i A;)), 
read 

Ao = , = 

A, = A^p- , A3/2 = ^ 

A _ n A — [(37r-128)/32+{256-187r)/34]7rV2 [b.b6) 

A2 - , A5/2 - 

A _ n A — (27/32-100/34+80/36)7r3/2 

As - , A7/2 - 24576R5 ■ 

We perform the subtraction proposed in (5.21) and we obtain the final result 

— r mim^R] + 



+ i6;^ + i2^J ^^("^^^^ 

/3^me 1 

16i? 1 - e^'^'"^-^ ■ ^ ■ ^ 

The functions Pn(a;) sue given by 



— 
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P3[X) 



P4:[X) 



dm 
dm 



m 



1-e 



27rm 



' 1\' 



m 



(5.65) 



they are all convergent integrals which can be easily numerically calculated. The 
finite part of the ground state energy given by (5.47) can be integrated by parts 
giving 



CSC 



oo 

-y 

m=l 

--/ 



dk k 



P 



n t 



dk k 



\nfo{ik) - 



0^ 
2kR 



(5.66) 



where the coefficient X„ , are given by (5.51), t = 1/(1 + {kR/mf)2 and 

hif^iik) = hi[CKm-p + ImK'^^p] + ln[/;i^^+^ + /„i^;^.J , (5.67) 

ln/o(2A;) = In [/^;^_^ + Joi^l^] . (5.68) 

Equations(5.64) and (5.66) are to be considered the main analytical result of 
this chapter devoted to the scalar field. Their sum gives the total renormalized 
vacuum energy. The sum will be performed in the numerical part of this chapter. 



5.5 Numerical evaluations 

In this section we show some graphics of E'^^'^ and of the complete renormalized 
vacuum energy £^l^ as a function of the radius of the string. We calculate also 
the asymptotic behavior of £^1 and E^^'^ for large and for small R. Since we want 
to study here only the dependence on R and on /3, we set me = 1 for all the 
calculations of this section. 

As a first step we rewrite in a form in which the dependence on the relevant 
parameters is more explicit: 

= ^ [P^9i{meR) + P'g2{meR)\ , (5.69) 



with 



gi{x) = [ -pi{x) + ^^P2{x) - \pz{x) + ^,Pa{x) - 77^ + ^ 



^^^^K-j g^.v-y 24^"^^ 96x 16(1 - 6^'^=-) 
giix) = ( -^^^2(a;) + ^P3(a;) + 1^ j . (5.70) 



55 



(5.71) 



the asymptotic behaviour of the p„(x) functions for x — > is found to be 

pi{x) ~ -1/24 +0{x) , 
1 1 

P2{x) - — + -In x- 0.0575 +C(a;) , 

4:X 2 

Psix) ~ ^ In X + 0.442 + 0{x) , 
3 

P4(x) ~ -lna; + 1.826 +0(a;) ; 

the corresponding behaviour of the functions 5'n(2^) is 

^i(x) ~ 0.0317+ C>(x) , 

g2{x) ~ 0.0417+ e)(x) . (5.72) 

The logarithmic contributions have cancelled. This was to expect also from the 
vanishing of the heat kernel coefficient A2. £^1 is proportional to R~'^ for — > 
and for an arbitrary (5. For it! —> 00 all the Pn{x) functions fall exponentially and 

so does £11 . 

The finite part £f^ is also proportional to in the limit i? — > 0, it can be 
seen substituting k ^ k/Rin the integrands of expression (5.66) 



CSC 

Cf - 



00 



27ri?2 ^ 

m=l 



^27ri?2 



dkk 



dk k 



n t 



lnfo{ik)\k^k/R - ^ 



(5.73) 



where the Jost functions are now independent of R. We name the first addend in 
(5.73) and second addend £jq, in the plots we will display them separately. 
For large R we found numerically £j'^ ~ R~^, which is in agreement with the 
heat kernel coefficient shown in (5.63), in fact the first non vanishing heat 
kernel coefficient after A2 determines the behaviour of the rcnormalized energy 
for R ^ 00. Below we show the plots of all the contributions to the renormalized 
ground state energy. The functions gl{x) and gl{x) are shown as well. Each 
contribution has been multiplied by R'^ so that all the curves take a finite value 
at i? = 0. We found necessary to sum up to 20 in the parameter m an to integrate 
up to 1000 in the k variable in order to obtain reliable plots. All the calculations 
where performed with computer programming relying on a precision of 34 digits. 



5.6 Discussion 

In this chapter we have carried out a complete calculation of the vacuum energy 
of a scalar field in the background of a magnetic string with delta function pro- 
file. The renormalized vacuum energy is given in terms of convergent integrals 
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X 



Figure 5.1: Scalar field. The functions gi{x) and g2{x) contributing to the 
asymptotic part of the ground state energy. 
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Figure 5.2: Scalar field. S^^,, E'L and E% multiplied by ■ for /5 = 2.2. 



(eq.(5.64) and (5.66)). A first remark can be made about the vanishing of the 
heat kernel coefficient A2. This coefficient, contributing to f^Jj'^^, is not zero for a 
generic background potential. The vanishing of A2 is also observed in a dielectric 



spherical shell with a squared profile |31|. Here, as in the preceding chapters 



it could be argued that singular profiles possess less ultraviolet divergences than 
smooth profiles. Such a statement is also confirmed by the heat kernel coefficients 
calculated in . 

The dependence of the sign of the energy on the radius R of the string and 
on the potential strength /3 is non trivial. The energy is negative only for large 
values of the potential strength, while for (5 smaller than one, the energy shows 
a maximum (Fig. 5.5). For large R the vacuum energy is negative in the back- 
ground of strong fluxes /5 > 1 and positive in the background of weak fluxes 
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Figure 5.3: Scalar field. The contributions m = 1, 2. ..5 to the finite part of the 
energy multiplied by ■ j3~'^, for j3 = 2.2. 

f3 < 1. This strong dependence on the parameter (3 was not observed in paper 
where an homogeneous field inside the fiux tube was investigated. In fact 
the most relevant result of our calculation is that the energy numerically shows 
a dependence on for large /3. 



. 8 p 




Figure 5.4: Scalar field. The three contributions to the renormalized vacuum 
energy multiplied by R^ ■ for /? = 0.4. The curve £^q shows the contribution 
to the finite part of the energy coming from the orbital momentum m = 0, while 
displays the contribution coming from the sum of all other m's. 
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Figure 5.5: Scalar field. The complete renormalized vacuum energy £^l^[R) 
multiplied by ■ P~^, for different values of strength of the potential. 
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Chapter 6 

Vacuum energy in the presence of 
a magnetic flux tube: Spinor case 

From the toy model analyzed in the last chapter, we pass here to the more 
interesting situation of a spinor field in the background of a magnetic string. The 
mathematical setup is analogous to that of the preceding chapter. The quantum 
field theoretical problem must be reformulated for the Fermi-Dirac statistic and 
for a quantum field with an additional degree of freedom. The Jost function will 
turn out to be simpler as that derived in the scalar situation. Furthermore the 
contribution to the vacuum energy coming from the orbital momentum m = 
will not need to be treated separately. 



6.1 Solution of the field equation 

The analysis of a spinor field in the background of a cylindrical magnetic field 
with an arbitrary profile has been performed in [39|. The field equation for a 
spinor field 

m = ( ''["I ) (6.1) 



, 92[r) 

in the background of a translationally invariant potential, with delta function 
profile is 



PO - me dr 

m+l-/3e(/?-r) 



m-(3e{R-r) 



r 




(6.2) 



The reader is referred to paper [39] for a derivation of this equation. Let us find 



the solution to (6.2) for one component of the spinor[|. The decoupled equation 

^We are not interested here in the complete set of solutions to eq.(6.2), the solution for g2{r) 
is sufficient to find the Jost function of the scattering problem. The decoupled equation for the 
component gi{r) is however indispensable and it will be used later. 
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for the component g2 is 



_ [m-PB(R-r)f ^ ^^^^ _ ^ ^ | ^^^^^ ^ ^ _ ^^ ^^ 



where k = y^pg — m^. The solution in the region r < R is 

gi{r) = Jmikr). (6.4) 

and in the region r > R 



1 



here f^'^{k) and /^*"*(A;) are the Jost function and its conjugate related to the 
scattering problem for the spinor field. Solutions (6.4) and (6.5) have the same 
form as those found in the scalar case for the inner and outer space. However at 
r = R the field has not the same form, as we discuss below, owing to the presence 
of the term j6{R — r) in the field equation (6.3). 

The ground state energy of the spinor field in the background of the magnetic 
string is 

= -y E 'Kt) . (6-6) 

n,a,a 

where the minus sign accounts for the change of the statistics, and the €(n^a,a) are 
the eigenvalues of the Hamiltonian 

H = -i-f^-f^ {d^i - ieAi{x)) + 7°me . (6.7) 

The degree of freedom a accounts for the two independent spin states. As in the 
scalar case we calculate the energy for a section of the string. The ground state 
energy density per unit length of the string in terms of the Jost function is given 

by 

S'P'^ = Cs / dk {k^ -miy-'dklnf^'''{ik) . (6.8) 

m=-oo -^"^^ 

The renormalization scheme is the same we introduced for the scalar case. The 
expansion of the ground state energy in powers of the mass and the definition 

of S^f^^ are the same as in (5.28), (5.29), apart from a factor —1 coming from 
the change of the statistics. The heat kernel coefficients will be of course not the 
same, we call them Bn- The normalization condition (5.30) remains unchanged. 
The ground state energy is split into the two parts 

^r"=^ E J dk[e-ml] — [In fZ'^itk)- In C^-^P^-{tk)] (6.9) 
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and 



^aT" = CsJ2 ^^[^' - ^e]'"^^ 1^ - > (6-10) 

where /^*^**'*"(iA;) is the uniform asymptotic expansion of the Jost function taken 
up to the third order in m. In (6.9) the analytical continuation s — > has already 
been performed, while in it will be performed after the subtraction of the 
divergent portion. (6.10) is a finite quantity for s — 0. 



6.2 Matching conditions and Jost function 



The matching conditions for the field on the surface of the string are not the same 
as in the scalar case. More exactly the condition for the first derivative of the 
field at r = is different from the one in the scalar case. The field is free inside 
the magnetic cylinder i.e. independent from /S, then from eq.(6.2) and (6.4) we 
have 

777 

kgiir) + drgiir) - -gi{r) = (6.11) 
and outside the cylinder 



kg^{r) + drg^{r)- 



Oi 



m — P 



0, 



(6.12) 



the continuity of the field on the boundary is required as in the scalar case: 



gi{r)\r=R = g2{r)\r=R, 9i{r)\r=R 

therefore, combining (6. 11), (6. 12) and (6.13) we find 

idrgiir))\r=R-idrg?ir))\r=R = -gUr). 



(6.13) 



(6.14) 



Finally, the matching conditions at r = i? read 



9?{r)\r^R 



92[r)\r=R 
{drgi{r))\r=R - {drg2{r))\r=R 



^giir)\r=R 



Inserting in (6.15) solutions (6.4) and (6.5) we get the system 



J' - i 



kR'^'^ 



for positive m, and the system 



J' - i 
7 ^ i 



kR'^-'^ 



(6.15) 



(6.16) 



(6.17) 
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for negative m. In (6.16) and (6.17) the argument (kR) of the Bessel and of the 
Hankel functions has been omitted for simphcity and the prime symbol indicates 
derivative with respect to this argument. Much in the same way as we did in the 
scalar case we find the Jost functions on the imaginary axis 

nr^itk) = t^kR[I^Km-^^, + Im+iKm-p], m>0 (6.18) 

/r""W - t-^kRilrnKrn-p+l+Im+lKrr^-^], m<0. (6.19) 

The two Jost functions are identical apart form the sign on the exponent of the 
imaginary factor. However formula (6.8) for the ground state energy contains 
the logarithm of the Jost function and the derivative with respect to k, then the 
imaginary factor i"^^ which is independent of k will not contribute to the energy. 

For the calculation of the asymptotic Jost function we found more convenient 
to use instead of m the expansion parameter 1/ given by 

m + 1/2 rmfor m = 0,1,2... (ro(W 
— m — 1/2 rmfor m = — 1,— 2,... 

with V — ■■■ in both cases. Then the Jost functions (6.18) and (6.19) become 

f:{ik) = + m>0 (6.21) 

which can be expanded for large positive m 

f-{ik) = + , m<0 (6.22) 

which can be expanded for large negative m. We recall the formulas for the 
asymptotic expansions of the Bessel I and K functions for u and k equally large, 
given by 

~ -^expj ^ a;"5/(n,a,t)| , (6.23) 

where x = 1/z/, the functions Si{n,a,t), SKin,a,t) arc the same as in the scalar 
case and a takes the values ±| for the Bessel I function and ±| ±/3 for the Bessel 
K function. From these formulas the logarithm of the asymptotic Jost function 
can be easily calculated up to the third order and we define 

]^fas-spin^ik^ = J2^J,n- , (6-25) 

j,n 

where t — + {kR/vY)^ and the nonzero coefficients are 

1^3,3 = - /5V12 , n,5 = -7/3V8 + /?V4 , (6.26) 
F3,7 = 5/3V8 . 
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6.3 The asymptotic and the finite part of the 
energy 

The asymptotic part of the energy can be written, using resuh (6.26), as 



p 



(6.27) 



we calculate the sum over v with the help of the Abel-Plana formula for half 
integer variables which can be found in section 3.3 of the chapter devoted to the 
spherical shell. The case m = (i.e. v = 1/2) does not need to be treated 
separately. We have only the two contributions 



pspin 



and 



nspm 
^as2 



47r V s V mi 



16R 



3,7 



,(-m2-2T(2-s))S„,,(mei?) . 



(6.28) 



(6.29) 



The functions S„j(x) are given in the appendix. They correspond to the functions 
Tjnj{x) found in paper H for a generic smooth background potential. The only 
pole term is contained in i^^fi" and the term proportional to me will be subtracted 
as well as the pole term, thus (6.28) cancels completely with the subtraction of 
^diT- "^^^ heat kernel coefficients which we have calculated up to the coefficient 
i?4, read 



Bo 


= 


) Bl/2 




= 87r/52 


, B3/2 


B2 


= 


, B5/2 


B3 


= 


, B7/2 
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_ 5/38 

1281ij6 







(3/3^+2/34)^3/2 
64_R3 

(135/3^-68/3-*+16/3g)7r3/2 
12288i?5 



(6.30) 



After the subtraction of i^X*" asymptotic part of the energy reads 



cspm 



~R2 



qi{mR) + 



■3^+6^ 



+ 

the functions g„(x) are 
qx[x) = 



+ 



127ri?2 67ri?2 



127ri?2 



q^irrieR) , (6.31) 



X 1 + e 



2ttu 



— 
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du 



dv 



du 



1 + e^'^'' 1/ 

V 

1 + e^^*" 

1/3 



7^ 



l + e 



27ri/ 



-)'- 



x^ 



(6.32) 



It is interesting to note how in both scalar and spinor cases we found the final 
expression for the asymptotic part of the energy to depend only on even powers 
of (3. This was to expect for physical reasons, in fact, inverting the direction of 
the magnetic flux the ground state energy should not change. 

The finite part of the energy S^j^^"^ can be hardly analytically simplified. We 
can only integrate expression (6.9) by parts to obtain a final form which is suitable 
for the numerical calculation: 



fspin 



" _i Jm, 



dk k 



3,7 



\nf^{ik)-j2yn,r 



p 



where 



Inf^iik) - In 



n,3 



(6.33) 



+ In 



(6.34) 



6.4 Numerical results 

We study numerically the behaviour of equations (6.33) and (6.31) and of their 
sum, which expresses the total renormalized energy, as functions of the radius of 
the string. The asymptotic behaviour of the two contribution is first found. The 
asymptotic part of the energy is rewritten in the form 



cspm 



where 



ei(,x) = (24gi(,x) - Aq2{x) + 7^3(0;) - qi{x)) 
e2{x) = 2 {q2{x) - q^ix)) . 

The asymptotic behaviour of the qn{x) functions for x — > is 

qi{x) - 1/48 +0{x) , 

q2{x) ~ ^ In a; + 0.635 +C>(x) , 



(6.35) 



(6.36) 
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qsix) ~ - In x + 1.135 +C»(a;) , 
q^{x) ~ ^lnx + 3.906 + 0{x) 

Zi 

and the corresponding behaviour of the e„(x) functions is 



(6.37) 



ei[x) 
62(2;) 



2 + Oix) , 
-1 + 0{x) ; 



(6.38) 



therefore behaves hke R ^ for small values of R. For R 



,-R 



00 " falls hke 



The finite part of the energy can be rewritten as 



pspin 



00 

R2 



dk k 



rrieR 



3,7 



In f^{ik)\k^k/R -^Yn,, 



(6.39) 



which is proportional to ^ for — > and to ^ for — > 00 in agreement with 
the heat kernel coefficient S5/2. We give the plots of the functions ei(x), 62(2;) of 
gspin^ £spm ^£ S^^, for small and for large values of the potential strength. 
The same remarks which we made in the scalar case about the numerical limits 
of summation and integration and about the precision employed in the computer 
graphics are valid here. 

ei{x) 



X 




Figure 6.1: Spinor field. The functions ei{x) and 62(2;) contributing to the 
asymptotic part of the ground state energy. 



6.5 Discussion 

In this chapter the vacuum energy of a spinor field in the presence of a semi- 
transparent magnetic was calculated, analytically in the first part of the chapter, 
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Figure 6.2: Spinor field. The curves of the asymptotic and of the finite part of 
the energy multiphed by R'^ ■ for /3 = 0.4. 



and numerically in the remaining part. As it is clear from the plots, the renor- 
malizcd vacuum energy changes with the radius, being not definitely positive or 
negative. 

With respect to the scalar case, we have an opposite behaviour in the region 
R < 1 (Fig. 6.5), where the energy shows a minimum for small values of /3, being 
positive for /3 > 1. However, when the radius become large, the vacuum energy 
shows the same behaviour in the scalar and in the spinor case: it is negative for 
large (3 and positive for small /3. In the limit /3 — > oo we found the energy to be 
proportional to /?^. This opens the interesting possibility that the vacuum energy 
could overwhelm the classical energy of the system when the fiux is sufficiently 
strong, in fact we have 

EroT - ^ + ^_g^ ' (6-40) 

^class ^vacuum 

where A, B and C are numbers and a is the fine structure constant. This pos- 
sibility is, however, not directly applicable to the model studied here. In fact, 
owing to the delta function the classical energy has not a finite value. 
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R 



. 02 L 



Figure 6.3: Spinor field. The curves of the asymptotic and of the finite part of 
the energy multiphed by ■ j3~^, for j3 = 2.2. 
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-0.02 L 

Figure 6.5: Spinor field. The complete renormalized vacuum energy S^p^^{R) 
multiplied by R'^ ■ /?~^, for different values of strength of the potential. 
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Chapter 7 
Conclusion 



We have studied the vacuum energy in the presence of curved semi-transparent 
boundaries represented by delta functions. The model and the results give rise 
to a number of discussions in the context of the existing literature. First, the 
sign of the Casimir energy in some of the plots provided in this work (Fig. 3.5, 
chapter 3, Fig. 6.5 chapter 6) is not in agreement with that found in calculations 
with similar models ||4^, ^ . Secondly, it seems arduous to extrapolate a general 
rule for the sign of the Casimir energy in the presence of a non-ideal boundary 
like a non-ideal sphere or a non-ideal cylinder. More exactly, in almost all of 
the situations investigated the sign of the energy is not fix, but it varies with 
the dimensions of the cavity. This cannot depend simply on the peculiarity 
of the potential profile chosen (which is singular and therefore unphysical) in 
fact, in papers [^, |7^, the vacuum energy of a scalar field in the presence of 
a spherical shell with Dirichlet boundary conditions and of a spinor field in a 
bag, was calculated and the sign was also found to depend on the radius. This 
result may express the interesting feature that the Casimir energy is in general 
a non monotonous function. It is also interesting the possibility that the energy 
could have a minimum in some situations, like it was found in chapter 3, for 
spherical shells with attractive potential or in chapter 6 for the magnetic string. 
This would imply that the Casimir force exerted on a hollow non-ideal object 
could vanish for some particular size of the object, and that the object could 
find some equilibrium in vacuum. However one should be careful in generalizing 
the results found here. We remeber that we have dealed for 3/4 of the work 
with the toy model of a scalar field, which does not have a counterpart in reality. 
Furthermore, in all the analyzed systems the classical energy has not been taken 
in to account. The classical contribution, which must undergo a renormalization, 
could compensate for the above mentioned effects. Only in the background of 
the magnetic string we found that the vacuum contribution could dominate the 
total energy of the system. However, to investigate extensively this aspect one 
should apply the method to a inhomogeneous magnetic string with non-singular 
potential, possibly taking also into account the anomalous magnetic moment. As 
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a matter of fact, the results found in chapter 5 are far from those found in 
where a homogeneous flux tube with square well profile is considered with much 
the same procedure employed here. In that work the vacuum energy was found 
to be negative in all the R-a.xes. This does not take place in our computations. 
However, in the concluding section of paper it was supposed that the presence 
of a singular background could considerably change the results. 

A common feature of all the investigated backgrounds is that the vacuum 
energy shows very weak UV-divergences. A crucial term in Sdiv, namely that given 
by the heat-kernel coefficient A2 was found to depend only on the third power 
of the coupling, while all the lower powers have cancelled. We had consequently 
a small number of pole terms in the calculations. The heat-kernel coefficients, 
which we calculated in the problem of the cylinder and of the magnetic string up 
to the coefficient A5/2, are confirmed by those calculated in other works [31, 35]. 

In the background of the magnetic string we found the coefficient A2 to be 
zero. Here the renormalized vacuum energy possesses no logarithmic behaviour. 
A logarithmic behaviour for small values of R was found in the background of a 
spherical shell (eq.(3.46)) and of a cylindrical shell (eq.(4.44)). 
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Appendix A 



Asymptotic expansion of the 
modified Bessel functions 



For the problems of the cylindrical shell and of the magnetic string, the expansions 
of the Bessel functions in negative powers of the parameters /, m and v (chapters 
4, 5 and 6, respectively) are obtained from the following formulas 

1 ^ 

Ii+a{kR) - ^=exp{ V r"5,(n,a,t)} , (A.l) 
Ki+a{kR) ~ J^exp{^r'^^;,(n,a,i)}, (A.2) 

n=— 1 

where a takes the values 1 and in the case of the cylinder and the values 
specified in sections 5.3 and 6.2 in the case of the string. The functions Si{n, a, t) 
and Sxin, a, t) are Deby polinomials. Up to the third order they are given by 



5*7(0, a, t) = -Int — -In 



2 2 \l-t 

Si(l,a,t) = -^{-3 + 12a^ + 12at + 5t'^) , 

Si{2, a, t) = -T^i^a^t + I2a^{-1 + 2f) + a{-26t + 30t^) + 3(1 - 6^^ + 5t^)] , 
48 

Si{3,a,t) = — (((25 - 104a^ + 16a^)f^)/3 + 16a(-7 + 4a2)t^ 
128 

-(531/5 - 224a2 + 16a^)t^ - (32a(-33 + 8a^)t^)/3 

-(-221 + 200a^)f - 240at^ - (1105t^)/9) ; (A.3) 

SK{-l,a,t) = -r^-lln(^i^^ , 
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SK{0,a,,t) = -lnt + -ln ^ 



2 2 \l-t 

+2 

Sk{2, a, t) = -j^i-^ci^t + I2a\-1 + 2t^) + a(-26i + 30i^) + 3(1 - 6^^ + 5i^)] , 
48 

5^(3, a, t) = — ^-((25 - 104a^ + 16a^)i^)/3 + 16a(-7 + 4a^)i^ 
128 

-(-531/5 + 224a^ - 16a^)t^ - (32a(-33 + 8a^)t^)/3 

-(221 - 200a2)t^ - 240at^ + (1105t^) /9 , (A.4) 

where t = (1 + {kR/Vf)~^ for the problem of the cyhndrical shell, t — {1 -\- 
{kR/m)'^)~2 for problem of a scalar field in the background of a string and t — 
(1 + {kR/vf)—2 for a spinor field in the background of a string. 
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Appendix B 

Calculation of the integrals 



In chapters 4, 5 and 6 the contributions to the asymptotic part of the energy 
coming from the Abel-Plana formula are calculated by means of the formulas 
given below, which were first derived in ||39|. The integration variable / is that 
for the cylinder. In the problem of the magnetic string it corresponds to m (scalar 
case) or to u (spinor case). The following formula 

2 2M-.^t^'_ m2-2^r(2-.)r(i±p)r(. + ^) 



dl / dk {k' -my-'dk — 



7" 2 (i?m)"-ir(j72) 

(B.l) 

has been used to obtain equations (4. 31), (5. 58), (5. 59) and (6.28). Formula 



r (I) (Rmy (1 + (^) 

has been applied to arrive at the results (4.32), (4.35), (5.62), (6.29). The func- 
tions Zn,j, A„j(x) and S„j(x) calculated for n < 3, j < 7 read 

4 r dl 



ZsAx) = +4/ dl{ -^—;\ Vl^-x^ 



4 

^3,5(3;) = +— 



4 r I ( f f p \' n' 1 



= dl-\[[ —-- - - v/p^ ; (B.3) 



15 \ VV 1 - e2^' / / / / 
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